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ABSTRACT 


The response of a constant- temperature hot-wire anemometer to 
sinusoidal and distorted sinusoidal acoustic oscillations is examined. 

The output of the anemometer is dependent upon the Reynolds number of 
the flow over the wire. The response is a measure of the interaction 
between the anemometer output and the acoustic pressure in the neighbor- 
hood of the wire. It is an open-loop prediction of the characteristics 
of actual closed-loop operation of a system. If the open-loop response 
is large enough, unstable closed-loop operation is predicted. 

The study was motivated by a need to investigate the stability 
limits of liquid-propellant rockets when perturbed by pressure 
oscillations. This is a very difficult process to study experimentally. 
The combustion process in a liquid-propellant rocket can be modeled as 
being Reynolds number dependent. The sinusoidal and distorted 
sinusoidal acoustic oscillations used for this study are the same as 
those characteristic of unstable rocket combustion. 

Both an analytical and an experimental study of the anemometer's 
response was performed. Qualitatively, the results of each are similar-- 
the response of the system to pure sinusoidal acoustic vibration of the 
fluid surrounding the wire is small, even when the magnitude of the 
acoustic pressure is quite large; but the response can be increased by 
as much as an order of magnitude with respect to the sinusoidal case 
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by the addition of distortion. The amplitude and phase of the distortion 
component, relative to the fundamental component, are the dominant factors 
in the increase in the response. 
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Chapter 1 


INTRODUCTION AND BACKGROUND 
INTRODUCTION 

The need exists for a more accurate understanding of the nonlinear 
instability limits of rocket engines. These limits are currently one of 
the least understood properties of combustion instability. The causes 
of instability and the factors determining the instability limits are 
of major importance to the design of high performance rocket engines. 

Setting 

Steady combustion in liquid propellant rockets is assumed to be 
vaporization limited -- a velocity sensitive or Reynolds number 
dependent process. Priem and Heidmann [1] have shown that quantitative 
agreement exists between theory and experiment for this mode. Vaporiza- 
tion limited combustion is also expected for unstable combustion and 
this has been qualitatively confirmed. 

Recently, Heidmann [2, 3] has shown analytically that distorted 
acoustic oscillations affect the open-loop response of vaporization 
limited combustion. The open- loop response is a prediction of the 
stability of the process during closed-loop operation; it indicates 
whether the acoustic gains of the system are sufficient to overcome the 
acoustic losses and cause instability. Response factors that were an 
order of magnitude greater than those obtained for pure sinusoidal 
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excitation were found for some conditions, and these response factors 
were high enough that unstable combustion was definitely predicted. 

The validity of Heidmann's results has not been tested experi- 
mentally due to the difficulties associated with making measurements 
in a burning rocket, the very short period between the onset of 
instability and the destruction of the rocket and the concomitant 
existence of many other uncontrolled variables. 

For the vaporization model used in references [2] and [3] it was 
assumed that the burning rate was proportional to the drop Reynolds 
number raised in a power. Such a model is characteristic of other 
processes such as those involving convective heat and mass transfer; 
hence, it appears that the same type of analysis could be performed for 
these processes. 

Scope 

In an effort to extend the significance of Heidmann's analysis 
and to confirm the results experimentally, another Reynolds number 
dependent process was selected and analyzed in the same manner as that 
done by Heidmann, and an experimental system was designed and used to 
test the results. The choice of the analogous process was primarily 
influenced by its amenability to experimental investigations. To 
facilitate comparison both the analytical and experimental studies, the 
imposed acoustic oscillations were the same type as those used by 
Heidmann — pure sinusoidal and distorted sinusoidal spinning transverse 
acoustic oscillations. Heidmann was particularly interested in the 
distorted sinusoidal oscillations since they are the most frequently 
encountered mode of unstable combustion [3]. Variations in the overall 
intensity and the proportion and phase angle of the distortion component 
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were considered. In addition, changes in the through-flow velocity were 
considered analytically but not experimentally. The open-loop results were 
then used to assess some of the conditions and mechanisms necessary for 
sustained closed-loop operation. In both Heidmann's analysis and this 
analysis, it was desired to see if the addition of distortion to sinusoidal 
acoustic perturbations of the Reynolds number dependent process involved 
would increase the response and drive the process to instability. 

BACKGROUND 

There has been considerable interest in studies of rocket combus- 
tion instability, expecially by NASA and the military. Combustion 
instability has been a problem to all modern rocket engine development 
programs from the German V-2's to the present [4]. At this time there 
is no way of reliably predicting whether or not a rocket design will 
be stable or unstable before actually building and testing it. Some- 
times during the operation of a developmental rocket the system will 
become unstable. Vibrations build up which, if allowed to continue, 
could destroy the rocket. These vibrations are easily detected by the 
loud screeching or screaming noise generated. 

Retrospect 

Up to the present time, the design procedure has been to create a 
rocket which is predicted to give the desired performance, then to test 
it to see if it will operate stably. If it does not, various modifica- 
tions are available which may be able to create stability. The primary 
methods are the insertion of baffles or acoustic liners into the combus- 
tion chamber [5]. These devices may damp out the generated vibrations 
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or prevent them from occurring at all; however, they may also have other 
effects on engine performance which may not be desirable. 

Many types of tests are used to check a new design. Scale model 
testing can be done; however, results from these tests are not always 
reliable. Sometimes the small model will be stable but the full-sized 
version will not be, and vice-versa. Static tests are also used in which 
the rocket engine is held stationary during operation while small explosions 
are set off inside the combustion chamber. In most cases the engine 
will remain stable up to a particular level of explosion, and self- 
sustained vibrations are triggered above that. These results again are 
not completely reliable. Sometimes low level explosions will trigger 
instability when more powerful ones will not. The type of explosion and 
its location also affect the results [3]. Hence, the only reliable way 
to determine stability is to actually flight test the full scale model. 

Since NASA sometimes has a human payload, the need exists to get as 
close to absolute reliability as possible. This becomes very expensive 
in both time and money; it would therefore be very desirable to 
develop some insight into why instabilities occur so they can be avoided. 

To try to better understand the causes of instabilities and the 
effects of various design parameters on them, many different types of 
studies have been initiated. These studies include the combustion 
process itself, the means and rate of fuel vaporization, the methods of 
injecting fuel and oxidizer into a liquid propellant system, the effect 
of the shape of solid fuel, the fuel/oxidizer ratio, the effect of 
nozzle configuration, the relative dimensions, and many others. Attempts 
have also been made to write the governing equations for the processes 
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inside a rocket; however, the system is so complex that no method has 
been satisfactory [6]. 

An approach to an understanding of the problem is suggested in 
some of Lord Rayleigh's work, which was originally published in 1878. 

He said that if an oscillatory heat source had a component of its rate 
in phase with an acoustic oscillation, the acoustic oscillation would 
be amplified [7]. This statement, coupled with the observation that 
certain definite patterns of vibration occur during instability, leads 
to a study of the interaction of the combustion process and the natural 
acoustic modes of vibration of the system. 

A liquid propellant rocket combustion chamber is basically a 
closed-ended cylinder through which a flow can pass. There are small 
holes in one end of the chamber for the injection of liquid fuel and 
oxidizer. At the opposite end of the chamber there is a nozzle 
through which the gases resulting from combustion pass out of the 
system. The nozzle is supersonic so that what goes on downstream of 
the throat section does not affect the chamber; however, there is a 
coupling between the nozzle and the chamber. Because of a rocket's close 
similarity to a closed-ended cylinder, this configuration has been used 
in this investigation. 

Inviscid Acoustic Fields of a Cylindrical Enclosure 

A knowledge of the characteristics of the solution of the inviscid 
wave equation for a cylindrical enclosure is necessary to an under- 
standing of both the theoretical and experimental analysis. The 
assumptions, boundary conditions, equations and details of the solution 
are given in Appendix A. The solution is a modification of the 
presentations given in references [8] and [9]. 
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The solution yields the frequencies of the natural modes of 
vibration of a closed-ended cylinder and determines the equations for 
the pressure, velocity, and density variations when excited at these 
frequencies. The equations show that there are three basic modes of 
vibration, each having an infinite number of permissible resonant 
frequencies. Also, one of these modes can occur in either of two ways — 
as a standing wave or as a spinning wave. Combination modes are also 
possible. 

The pressure, velocity and density can each be expressed as the 
sum of a time-mean and a time-dependent component 

P(r,e,z,t) = P(r,e,z) + P'(r,e,z,t) 

V(r,e,z,t) = V(r,e,z) + V'(r,e,z,t) (A. 7) 

p(r,s,z,t) = pT(r,e,z) + p'(r,0,z,t) 


where the bar superscript denotes the time-mean value and the prime 
denotes the deviation from the mean value. From Appendix A the 
equations for the time-dependent components of pressure, velocity, and 
density are 


P' = p* 
mnn z mnn z 


J (a irr/R) 

‘ ' cos(n,Trz/L)cos(me)sin(iu „ t + B ) 

z mnn z mnn z 


m run ' 


(A. 44a) 


P mnn “mn^j' (a m _irr/R) 


z - m mn cos(n_rrz/L)cos(me)cos(a) mnn t+e mnn ) (A. 44b) 


V ' j r 

r „ p to__„ R J m (a „ir) 
mnn z «> mnn z nr mn 


mnn z mnn z 


¥ i 


mnn_'" 0 m (a 7rr/R) 


P m 

z ~nT"mn 
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The first pure mode of vibration is a longitudinal type for which 
m and n are zero. The resonant frequencies of this mode are whole 
multiples of the fundamental frequency of this mode. The pressure and 
velocity distributions for the first three resonant frequencies at quarter 
cycle intervals of time are shown in figure 1; The drawings above the 
pressure distributions combine the pressure and velocity into one 
picture. The denser the dots, the greater the pressure and density of 
the gas. The arrows indicate the net motion of the molecules in both 
direction and magnitude. Both the pressure and velocity vary sinusoidally 
along the length of the chamber, but are ninety degrees out of phase with 
each other in time and space. There is no variation in pressure or 
velocity either radially or tangentially. 

The next two modes of vibration are transverse types where the 
variations in pressure and velocity can occur tangentially and/or 
radially but not longitudinally. The resonant frequencies of these 
modes are not whole multiples of a fundamental. The frequencies are in 
part determined by the zeros of a Bessel function, which is not a true 
periodic function. 

The first transverse type is a pure radial mode; there is no 
tangential variation in pressure or velocity. For this mode m and n z 
are zero. The pressure distributions at wt = 0 and it for the first three 
radial modes are shown in figure 2. The radial variation in pressure 
is a first order Bessel function; the variation with time at a particular 
location is sinusoidal. Figure 3 includes the velocity distribution 



Figure 1 . 
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and drawings which pictorially show how the pressure and velocity vary. 
The maximum velocities occur at pressure nodes. The overall velocity 
effect is movement of the gas molecules in and out along radial lines. 

The pressure and velocity are a quarter cycle out of phase with each 
other in time. 

The third mode of vibration is the tangential one which has both 
radial and tangential variations in velocity and pressure. For this mode 
n and n z are zero. It can occur in two forms -- standing and spinning. 
The standing tangential is illustrated in figure 4. The pressure distri- 
bution is a Bessel function -- the first standing tangential having a 
first order Bessel function variation, the second having a second order 
Bessel function variation, etc. The pressure in the first varies from 
high on one side of the chamber to low on the other side, then reverses 
itself. This causes the gas in the chamber to "slosh" back and forth 
across the chamber. When the pressure distribution has reached a 
maximum, the velocity distribution is zero, and vice versa. The pressure 
and velocity are ninety degrees out of phase with each other. 

In the second standing tangential the "sloshing" is between 
quadrants rather than from one side to the other. As the mode number 
increases, the number of pressure nodal lines increases; hence, the 
number of segments into which the chamber is divided becomes greater. 

The spinning tangential mode has the same resonant frequencies, 
the same Bessel function distributions in pressure, and the same 
velocity distributions as the standing one; but their orientation 
changes (spins) around the chamber at the resonant frequency. The 
equations have a slightly different form, namely 
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Figure 4. Pressure Distributions and Schematic Velocity Distributions for the First Three 
Standing Modes of Vibration (PN = Pressure Node) 
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There is an interchangabil ity in the effect of time and angular position. 
The standing and spinning tangential modes are compared in figure 5. 

There are no longer any pressure nodal lines -- there is a nodal point 
in the center. At any location in the chamber the maximum pressure 
and maximum tangential velocity occur at the same time -- they are in 
phase. But the pressure and radial velocity are out of phase by ninety 
degrees . 

What particular mode or modes of vibration will develop in a 
chamber, and how intense they will be, depends on the nature of the 
driving mechanism. The driving mechanism needs to be located in the 
vicinity of a pressure antinode for the particular mode being driven. 
Referring to figure 1, the driver should be placed near either end to 
create the first longitudinal mode. To create the second longitudinal 
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mode, the driver could be placed near either end or the middle. To 
confirm that the location of the driving mechanism affects which mode 
will develop in a rocket, experiments have been performed in which 
various modes of instability were generated by changing the location of 
propellant injection [10]. Also, examinations of rockets which have 
developed an unstable pattern of combustion during actual operation 
indicate that strong spinning transverse modes of acoustic oscillations 
are dominant. 

An effective location for pressure and velocity probes can also 
be determined by the pressure and velocity distributions. For example, 
a pressure probe at the center of the chamber would not detect any 
variations due to the first longitudinal mode since that is the location 
of a pressure node; however, a velocity probe at the center would 
detect maximum velocity variations for the first longitudinal mode. 


Response Factor 


As was mentioned earlier. Lord Rayeigh suggested that the 


interaction between an oscillatory heat source (e.g., a combustion 
process) and acoustic oscillations could amplify the oscillations [7]. 


To measure this coupling between the combustion process in a rocket 


engine and the pressure variations associated with the acoustic 


vibrations, nonlinear, in-phase response factor has been developed 
[11] and is defined as 




( 1 . 1 ) 
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where W = perturbation in the normalized burning rate 

P' = perturbation in the pressure 
The notation being used here is of particular significance. As 
was done earlier, the conventions of acoustics are being used in that 
the value of a property is expressed as the sum of a time-mean value 
and a perturbation about the mean, i.e., 

& * ? + e' (1.2) 

where the bar over the property indicates the time-mean component and 
the prime designates the time-dependent component. Further, the time- 
dependent component can be expressed in terms of a non-dimensional 
component b 1 , 


B 1 = (1 >3) 

Using this notation, b 1 represents a percentage, or fractional, deviation 
about the mean value. Substituting equation (1.3) into equation (1.2) 
and regrouping, p becomes 

6 = B (1 + l') (1.4) 

It is the perturbation about the mean value s' which is of interest, not 
the mean value itself. It has been shown that unstable combustion 
for sinusoidal oscillations will occur when the in-phase response 
factor is greater than some value between 0.8 and 1.0; that is, the 
oscillations will grow with time [3]. A nonlinear out-of-phase response 
factor, I has also been developed, but its role is not entirely clear 
at this time. 
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Studies have shown that although vibrations in a combustion 
chamber may start out as sinusoidal, they change with time and develop 
into a distorted wave form [12]. Realizing this, Heidmann of NASA's 
Lewis Research Center analytically calculated response factors for 
distorted sinusoidal wave forms [3]. He used as a model of the combus- 
tion process 



where w = burning rate (energy release per unit time and volume) 

Ci = constant of proportional ity 
Re = propellant drop Reynolds number 
p = gas density 

D = diameter of the propellant drops 

U = magnitude of the relative velocity of the drops compared 
to the gas 
v - gas viscosity 

m * Reynolds number exponent which characterizes the 
particular type of burning 

This Reynolds number dependency is characteristic of many combustion 
mechanisms and in particular is appropriate for the drop vaporization 
process in a liquid propellant rocket [13]. Assuming that burning 
takes place immediately upon vaporization of the propellants, and that 
m = 0.5 is typical of the drop vaporization process, the burning rate 
for vaporization limited combustion, using equation (1,5), is 


pDU 

VO - c, 

1 u 


0.5 


( 1 . 6 ) 
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It is this expression for the burning rate which Heidmann used in his 
analysis. From it he created a normalized burning rate W such that 


W 


M? 


IV 


0.7) 


where to = burning rate with oscillations 

to o = burning rate without oscillations 

To expand W, an expression is needed for the relative velocity of the 
drops in terms of steady and oscillatory components. It was assumed 
that the droplets have only a steady axial component of velocity. Then 
the magnitude of the relative velocity of a drop compared to the local 
gas velocity is given by 


U 2 = | v 2 t 2 + v; 2 + v; 2 (1.8) 

where u £ = axial velocity of a drop 
v 2 = axial velocity of the gas 
Vg = tangential velocity of the gas 
v^, = radial velocity of the gas 

The axial velocity difference is defined as a V 


aV = |u £ - v z | (1.9) 

and is assumed constant. The velocity components v^ and v^. represent 
the acoustic perturbations of the gas; the perturbations occur only in 
the transverse directions. They are combined into one term v£ defined by 


v! 2 = v 1 2 + v 1 2 
t e r 


( 1 . 10 ) 
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From equations (1.8), (1.9) and (1.10), 



f V tl 

2 ' 

= AV 

1 + ki 



l AV J 



0.5 


( 1 . 11 ) 


Writing p in terms of its time-mean and time-dependent components, using 
equation (1.11) to express the velocity, and equation (1.5) for the 
burning rate* one obtains 



Simplifying, this becomes 


+ pT 

i + 

i| 2 ' 


. 

AV 1 . 


m /2 


( 1 . 12 ) 


The perturbations in density and velocity can be expressed as Fourier 
series as follows: 


p' = r p cos(na>t - 4> ) (1-13) 

n=l n n 

ca 

v t = z c » v n cos ( nu)t - * n * 0 n > 0-14) 

n=l 

where the <t> n ' s specify the phase relationships between the harmonic 
components, and the e n 's specify the phase relationships between the 
density and velocity. As was discussed in Appendix A, for a spinning 
situation e n = 0. 
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Heidmann devised these expressions to reflect the observation 
that the oscillations which occur during unstable combustion seem to 
be predominantly composed of a spinning tangential oscillation at a 
fundamental frequency distorted by a spinning tangential oscillation 
at the second harmonic frequency, plus other lesser contributions at 
higher harmonics. The expressions are modifications of those obtained 
from the solution of the wave equation, equations (A. 48). (The dependence 
on position does not appear in the equations since the calculations are 
carried out for a representative location. Also, the density has been 
normalized.) The major difference between Heidmann's equations and 
equations (A. 48) lies in the characterization of the frequencies. 

Heidmann characterized them as being whole multiples of the fundamental 
frequency since this was what was observed in actual rockets; however, 
the spinning tangential resonant frequencies obtained from the wave 
equation are not whole multiples of the fundamental due to their 
dependence on a Bessel function. This point will be discussed in more 
detail later. 

Substituting p' and v^. into equation (1.12), expanding and 
rearranging terms, an expression for the normalized burning rate can be 
obtained in the form 


W ■ W (1 + W') (1.15) 


W’ can be written as the harmonic series. 


00 CO 

w' s Z z a cos (nut + a ) (1.16) 

n—1 r=l n,r n,r 


where the A Vs are the harmonic coefficients for burning rate perturbations 
n s r 
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and the a 's are the harmonic phase angles for burning rate perturbations. 
n»r 

The pressure perturbations can also be expressed as a series, 

oo 

P' = z p n cos(nwt - <f> n ) (1.17) 

n=l 

W 1 and P* can then be substituted into equation (1.1) for and 
integrated to obtain the response. Instead of carrying out these rather 
involved calculations, numerical methods can also be used to obtain the 
response. (Heidmann used both techniques and found good agreement 
between them. ) 

The last part of this analysis is only outlined here since the 
actual expressions involved are very long and complex, and it is not 
the purpose here to give a detailed explanation of Heidmann 's work but 
rather to give an outline of his procedures and results so that they 
can be compared to the procedures used for this research. 

Heidmann calculated for various combinations of harmonic 
content, phase angles, relative drop velocity, and other variables. He 
found that the system is predicted to be relatively insensitive to 
sinusoidal perturbations; however, adding varying amounts of distortion 
could increase R^ by an order of magnitude. Unstable combustion was 
definitely predicted for some conditions. The relative axial drop 
velocity was also shown to affect the stability. 

Experimental confirmation of these results was desired but could 
not be obtained due to the problems associated with making measurements 
in a burning rocket, the very short period between the onset of 
instability and the destruction of the rocket and the concomitant 
existence of many other uncontrolled variables. 
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PURPOSE AND SCOPE 

The purpose of this work is to select a Reynolds number dependent 
process which is amenable to experimental investigation, to perform a 
Heidmann-type analysis of it to determine if it is qualitatively the 
same and to perform an experimental study for the same conditions examined 
theoretically — insofar as is possible within experimental limitations. 
The results of the above studies will then be compared and used to 
determine if the addition of distortion to sinusoidal acoustic perturba- 
tions of the process involved will increase the response and possibly 
drive it to instability. These open-loop results will then be used to 
explain the mechanisms necessary for sustained closed-loop operation. 



Chapter 2 


ANALYTICAL INVESTIGATIONS 
INTRODUCTION 

To extend the significance of Heidmann's analysis and to have 
a basis of comparison for the experimental work, another Reynolds 
number dependent process was selected and analyzed in the same manner 
as that by Heidmann. To facilitate comparison, the same type acoustic 
oscillations were imposed — pure sinusoidal and distorted sinusoidal 
spinning transverse acoustic oscillations. Heidmann was particularly 
interested in the distorted sinusoidal oscillations since they are 
the most frequently encountered mode of unstable combustion [3]. The 
choice of the analogous process was primarily influenced by its 
amenability to experimental investigation. In both Heidmann's analysis 
and this analysis, it was desired to see if the addition of distortion 
to sinusoidal acoustic perturbations of the Reynolds number dependent 
process involved would increase the response and drive the process to 
instability. Variations in the through-flow velocity and the amplitude 
and phase angle of the distortion component were considered. The 
analytical investigations presented are based on reference [14]. 

THE ANALOG 

l 

A qualitative description of the drop vaporization process will 
provide a basis for choosing a suitable analog. In the combustion 
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process a liquid droplet exists at a temperature close to its saturation 
temperature corresponding to the partial pressure of its vapor. The gas 
mixture, which surrounds the droplet, moves relative to it and has a 
temperature which is often significantly higher than the droplet's. 

Much of the energy transfer from the gas mixture to the droplet is by 
forced convection and, under steady-state conditions, results in 
vaporization of the droplet liquid. 

Assuming the drop vaporization process to be convection dominated, 
a reasonable analog of this process is that of forced convection heat 
transfer from a hot wire to a cool gas. Since the constant-temperature 
hot-wire anemometer provides an excellent means for experimentally 
determining the instantaneous energy transfer rate from a hot wire, it 
is used as the analog process. 

To provide a controlled unsteady environment for the hot wire, a 
reasonant sound field with forced mean gas motion is used. Before 
giving a detailed description of the wire's environment, the physics of 
a constant-temperature hot wire in an unsteady gas flow will be developed. 

The physical system is shown in figure 6. It is assumed that 
there exists a steady mean flow w, that unsteady perturbations in the 
flow may exist in any direction, and that the wire is oriented such 
that its axis is normal to the mean-flow direction. For the cartesian 
coordinate system shown in figure 6, the z-axis is in the mean-flow 
direction and the y-axis is parallel to the cylinder. 

Energy is electrically supplied to the wire such that its 
temperature is constant; the constant- temperature anemometer provides 
the needed energy and controls. A natural limitation for constant 




Figure 6. Schematic Diagram of the Hot Wire and Its Environment 
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temperature operation is that the gas temperature must not exceed the 
normal operating temperature of the wire. 

Taking the hot wire to be a closed system, the first law of 
thermodynamics is 


d£ = dE . dW 
dt dt dt 


( 2 . 1 ) 


For constant temperature operation the total energy of the wire is 
constant and equation (2.1) becomes 


dQ _ dW 
dt “ dt 


( 2 . 2 ) 


The rate of electrical work done by the system is 


3 ?=-'% < 2 ' 3 > 

Three modes of energy transfer as heat exist; however, radiation 
is small enough to be neglected, conduction and convection are lumped 
together, and Newton's law of cooling is used to express the heat transfer 
rate. Thus 


di = . h A [T - T] 
dt w L w J 


Introducing the Nusselt number. 


hD. 

Nu f = £ 


w 


equations (2.2), (2.3) and (2.4) reduce to 


I2R W = * k f L w N “f£ T w - 


(2.4) 


(2.5) 


( 2 . 6 ) 
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Certain auxiliary relations are now introduced to permit equation 
(2.6) to be expressed in terms of the properties normally measured and/or 
known. 

The electrical resistance of the wire is a linear function of 
temperature in the range being considered; therefore, its operating 
resistance can be expressed as a function of its "cold" resistance 
R ra and the temperature difference between the atmosphere and the wire 
when it is in operation 

«*= R . + cvo < 2 - 7 > 

Since in the Disa anemometer system the potential drop across the 
wire is not directly measured, it is necessary to express the current in 
equation (2.6) in terms of the Disa bridge properties. For the bridge 
network shown in figure 7, the hot-wire current is 


I = 


R + R . , + R + R , 
s cbl p w 


( 2 . 8 ) 


There are a number of possible correlations that could be used 
to express the Nusselt number Nu^ in terms of the flow field variables. 
King [15], in 1914, was the first to develop such a correlation. Recently 
Andrews, Bradley and Hundy [16] brought together the many different 
correlations which have been proposed by various investigators since then 
and compared them to each other and to some extensive work of their own. 
They found that the most reliable correlation for low Reynolds number 
flows (0.02 < Re^ < 44) is that of Col 1 i s and Williams [17], namely 


Nu f “ [A + B Re f N ][T f /Tj 0 ' 17 (2.9) 


where A, B and N are determined by actual calibration of the hot wire 
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used and Tf is the film temperature defined as = (T w + T ro )/2. (The 
calibration procedure is described in Chapter 3.) 

Since the experiments used atmospheric air at room temperature 
and pressure, the viscosity and thermal conductivity data are taken from 
Hilsenrath [18] and the density is computed from the ideal gas relation, 
namely 


p - 


R a T 


( 2 . 10 ) 


The Hilsenrath data in the temperature range of interest (273.15°K < T < 
473.15°K), can be accurately approximated by 


and 


k - k 0 [T/T 0 ] 


“ * ‘•owy 


0.877 


0.757 


( 2 . 11 ) 


(2.12) 


The relationship between the hot-wire bridge voltage and the flow 
properties is obtained by combining equations (2.6), (2.8), (2.9), (2.10), 
(2.11) and (2.12). This gives 


E 2 = 


't R S t R cbl + R p * V 2 


O L w k 0 ] [a + B Re f ^ N (V 2 ) N/2 ] 


0.17rj 1 0.877 


TT T 0.l7r T *1 
if l f 

It J ItJ 


[T - T] 
L w J 


(2.13) 


where 


Re 


P Vj D, 


w 


f n r- 1 - /-r- -.0-757 

1 R a T f p o^W 


(2.14) 
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and 


V 2 = [(u’) 2 + (Kv 1 ) 2 + (w + w 1 ) 2 ]/V i 2 (2.15) 

K represents the sensitivity of the hot wire to flow parallel to its 
axis and is approximately equal to zero. (The wire is primarily sensitive 
to the components of the flow normal to its axis; it is a two- 
dimensional sensor. Also, the wire cannot distinguish direction. A 

~y 

velocity of V = Vi generates the same response as a velocity of V = Vk. 
This means the wire cannot distinguish between positive and negative 
flows either. A velocity of V = V Q sin at i will appear to the wire 

as the rectified speed, V = |V Q sin cut|.) 

Since the empirical relationship for Nu^, equation (2.9), was 
obtained for steady-state, steady-flow conditions, its applicability to 
unsteady flow may be questionable; however, it is customary to assume 
that it holds for unsteady flows if instantaneous values are substituted 

for the fluid properties [19]. This is done here and it constitutes a 

major assumption. 

Noting that the leading coefficient in equation (2.13) is constant 
during normal operation and that it depends on the wire temperature, it 
can be rewritten in a simpler form by defining dimensional forms of A 
and B, namely 

A = C(T W )A (2.16a) 

B = C(T w )B (2.16b) 


where 
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C(Tj = ^k, k n 
' w wo 


K * R chl * *0 * U ! 


w 


With this nomenclature the final form of equation (2.13) is 


■[ 


A + B Re. 


. _ T -0.17 rT -.0.977 ">0.5 

<* 2 > N/ 2 ][r] [r] [ V T i; 


(2.17) 


E is a function of the pressure, velocity and temperature; to determine 
E the wire must be calibrated and the flow field properties must be 
known . 

During the steady-state, steady-flow calibration of the wire all 
quantities in equation (2.17) are constant except the velocity; hence, 
equation (2.17) reduces to 


where 


E * [a + B V N ]°’ 5 

4T.-lO.17fT 10,877 

A = A lr] LrJ £t w - t] 


(2.18) 


(2.19a) 


and 


N[!f]°- 17 [!i] 0 - 877 


8 = B R 


[T w - T] 


(2.19b) 


The calibration determines the constants, A, B and N, from which A 

A 

and B and then A and B can be determined. 


FLOW FIELD PROPERTIES 


The wire's environment and orientation are schematically shown 
in figure 8. (The wire is no longer oriented parallel to the y-axis 




Figure 8. Schematic Diagram of the Hot-Wire In its Cylindrical Enclosure (One End Plate Is Not Shown) 
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as it was in figure 6.) The wire is surrounded by a cylindrical enclosure, 
with inside diameter D and length L, which is acoustically closed but 
permits a uniform forced flow in the axial direction, (This can be 
accomplished experimentally by the use of perforated end-plates.) The 
wire is located in a transverse plane. 

With proper excitation of the enclosure, acoustic velocity 
oscillations in any of the three vector directions can be obtained; 
but, as was described in Chapter 1, Heidmann [3] restricted his analysis 
to transverse-velocity oscillations "because of the greater interest in 
transverse-mode instability in rocket combustors." Considering only 
transverse oscillations to exist, the velocity at any point in the chamber 
can be expressed as the sum of a constant axial velocity 7 Z and an 
oscillating transverse velocity v^., i.e., 

V 2 = v z 2 + v.j. 2 (2.20) 

where 

v’ 2 = v r 2 + v e 2 (1 ‘ 9) 

The wire is oriented such that it makes an angle a with the radius vector 
(figure 8). As has been stated previously, the wire is not sensitive 

to flow parallel to its axis, therefore the wire cooling velocity V c , 

expressed in terms of the enclosure's cylindrical system, is 

V 2 = v, 2 + (v' sin a + v ' cos a) 2 
c z r y 

Nondimen sional izing this with V x 2 , gives 

V c 2 - [v z 2 + (v; sin a + Vq cos a) 2 ]/!^ 2 


( 2 . 21 ) 
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For the coordinate system shown in figure 8, V c 2 is used instead of V 2 
in equation (2.17), the equation for E. 

The equations for the acoustic pressure and particle velocity 
are those obtained in Appendix A, equations (A. 48). The acoustic 
temperature is related to the acoustic pressure. Since the process is 
assumed to be isentropic. 


p ( 1 -y) t y = b (2.22) 

Expressing P and T in terms of their time-mean and fluctuating components, 
equation (2.22) becomes 

(P + P')^" y) (T + T’) y = b (2.23) 

Expanding the terms on the left-hand side of equation (2.23) by use of 
a binomial expansion, using equation (2.22) for b and eliminating higher 
order terms yields 


Y — + (1-y) — = 0 (2.24) 

7 P 

In terms of P' and T' , equation (2.24) is 

T' - i 1 —) P' (2.25) 

The equations for the oscillations in pressure, velocity and 
temperature are 


P' - 
mnn 


z 


J n/ 0 W nr >'R) 

= L n — i - 7 s — cos(n _n-z/L)sin(me-w t-n n ) 

mnn z z mnn z mnn z 


(A. 48a) 
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(a. 48b) 
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(A. 48c) 


(2.26) 


where T = T and P = P . 


NONLINEAR RESPONSE FACTOR 


Exposing a constant-temperature hot-wire to a mean flow perturbed 
by acoustic oscillations results in an unsteady hot-wire voltage E, 
Expressed in terms of its mean and fluctuating components, E is 

E = E + E‘ (2.27) 

or 

E = E[1 + E 1 ] (2.28) 

with 

E* = E 7E. (2.29) 

A nonlinear, in-phase response factor can be formed for this 
process just as was done for the combustion process, i.e.. 



(2.30) 
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It is a measure of the portion of E' which is in phase with the acoustic 
pressure P\ 

To evaluate E ’ (cut ) must be related to the sound field given 

by equations (A. 48a), (A. 48b), (A.48C) and (2.26) as follows: 

©Certain constants of the anemometer and the environment must 

be determined. These include R s » , R p P w , T^, k Q and y Q . 

©A hot wire is selected and its operating temperature determined. 

This establishes T , T , L and D . 

WWW w 

©The hot wire is calibrated under steady-state, steady-flow 
conditions. From the calibration, equation (2.18) is 
obtained. From this, the constants A and B can be calculated. 
©The hot wire's environment is specified. This includes 7 z> 

P' and the wire's location in the chamber. 

©All the above values are substituted into equation (2.17), 
thereby producing £(wt), (When the above quantities have 
been specified, the only variable remaining is wt.) 

©The time-mean value of E(tot) , E, is obtained from 

2tt 

E = ^ E((ot)d(<ot) (2.31) 

/ 

o 

Actual mathematical integration of E (cot ) is difficult", 
therefore, the integration is performed numerically. 

• The desired E'(wt) is given by 

E* (tot) = [E(«t)/E] - 1 (2.32) 

©With E' given by equation (2.32) and P' given by equation 
(A. 48a) nondimensionalized over P^, is obtained from (2.30). 



37 

As before, the integration of E{wt) is done numerically. 

The nature of the nonlinear response factor is best 
illustrated by numerical examples. 

NUMERICAL EXAMPLES 

The numerical examples to be presented are based on the following 
hot-wire and environmental properties: 

• Wire and Bridge Properties 

L = 1.02 mm 
w 

D w = 5 ym 
a = 0.0126 ohm/K° 

Rp = 0.40 ohm 
R , , = 0.30 ohm 

CD i 

R„ = 100 ohms 
s 

R - R =1.40 ohms 

w 00 

R = 3.59 

ot> 

A = 0.521 
B = 0.673 
N = 0.444 

• Basic Air Properties 

T q = 273.1 5°K 

k Q = 2.413 x 10~ 2 kg-m/(sec 3 -K°) 

y Q = 1.716 x 10‘ 5 kg/(sec-m) 

R a = 287.1 m 2 /(sec 2 -K°) 

y = 1.4 
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©Wire Environmental Properties 

T = 298. 3°K 

co 

= 735.0 mm Hg 
L = 31 .30 cm 
D = 20.49 cm 


With these properties specified, the main variables are the hot- 
wire Mach number M = v /c , the acoustic oscillations and the location 

z 00 

and orientation of the wire. 

The equations used to express the acoustic oscillations have 
to be compatible with those used by Heidmann and those obtained from 
the solution of the inviscid wave equation. As was mentioned in 
Chapter 1, Heidmann assumed the frequencies involved where all whole 
multiples of the fundamental frequency since that was what has been 
observed in actual rockets. He regarded the acoustic oscillations as 
the sum of a fundamental (n = 1) plus distortion components (n = 2,3,...). 
In his calculations he used only one distortion component since, in 
actual rockets, the contribution from the higher harmonics is very 
smal 1 . 

From equation (A. 38) 



“Win, c 
z “ 

2n 2 


ot 2 n 2~\ ® ^ 

(- x ■) + <-f> 


(A. 38) 


which expresses the resonant frequencies as derived from the inviscid 
wave equation, it is clear that the resonant frequencies of pure 
spinning tangential waves are not whole multiples of the fundamental; 
however, by proper choice of the L/D ratio of the chamber, a combination 
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mode can be forced to have twice the frequency of the pure mode. This 
was done here. A more complete discussion of this is given in Chapter 3. 

With L and D as given, the frequency of the second-spinning, second- 
longitudinal combination mode (hereafter referred to as the (202) mode) 
is twice the frequency of the first-spinning pure tangential mode (here- 
after referred to as the (100) mode). The (202) mode has a time-dependent 
component of longitudinal velocity v z '; however, a longitudinal velocity 
node exists for this mode at the central cross-section of the chamber. 
Because of this, the spatial behavior of the (202) mode, at the center 
of the chamber, is the same as that of a second spinning pure tangential 
mode (hereafter referred to as the (200) mode). Hence, the longitudinal 
location of the wire is chosen as z = L/2. 

To have agreement with Heidmann in the form of the velocity 
equations, the wire is chosen to be oriented such that a = 0 (the wire's 
axis is parallel to the radius vector) and e = 0. In this way only 
the tangential component of velocity affects the wire and the form of 
the velocity equation agrees with that of Heidmann. The radial position 
of the wire is taken as r = 0.95R. 

When the above quantities are substituted into the equations for 
the acoustic pressure and particle velocity, they have the same form as 
Heidmann 1 s equations for the acoustic pressure and particle velocity. 
Letting P{ be the normalized acoustic pressure of the (100) mode and 
be the normalized acoustic pressure of the (202) mode, and with 
w 2 02 = 2w 10 o = 2u>, the fundamental and distortion components of the 
pressure are 


p i = p ioo/ p „ = Pi cos (<ot - <h) 


(2.33a) 
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and 


P 2 — ^202^00 ~ P2 COS(2oit - <}> 2 ) 

where 

PlOO 0i(0.95 a 10 it) 

Pl P*, Jl"( a l0 71 7 

*1 s - (nioo + 90°) 

-A 

P202 02(0.95 020 TT ) 

P2 P*, 0 2 (a 2 o Tr) 

and 

<J>2 = "(^202 + 90 °) 

Similarly, the acoustic particle velocities are 

v + ‘ * v ' = c Vi cos(cot - 4 >! ) 

tj 8 100 


and 


V = v ' = c v 2 cos(2u)t - * 2 ) 

t 2 e 202 


where 


C Vi = 


100 


.1. (n qc; 


and 


2 P 


c v 2 = 

00 4 


202 


J 2 {0. 95 a 2 o tt) 


0.95 p (0202 R ^ 2 V^T 


(2,33b) 

(2,34a) 

(2.34b) 

( 2 . 34c ) 

(2.35a) 

(2.35b) 

(2.36a) 


(2.36b) 
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Combining equations (2.33) and equations (2.35), gives 

2 

P' = z p'. cos (nut - $ ) (2.37a) 

n=l n 

and 

2 

v! = z c v„ cos(nwt - ) (2.37b) 

t , ® n n 

n=l 

Equation (2.37a) is the same as Heidmann's equation (1.17) with the 
series terminated at 2, as was done in his analysis. Equation (2.37b) 
is the same as Heidmann's equation (1.14) with the series again terminated 
at 2 and with e n = 0, which again is what Heidmann used. (The 0 n ‘s were 
set equal to zero by Heidmann to model a spinning situation.) 

In carrying out the actual calculations of R^, the wire Mach 
number (M = v /c ), the intensity of p x , the ratio of p 2 /Pi and the phase 

Z 00 

angle between p 2 and p x can be varied. What follows is a presentation 
of the effect of these variables on R^. 

RESULTS 

The following is a graphical presentation and discussion of the 
results of the calculations of the response factors for various values 
of M, Pj, p 2 and <f> 2 

Influence of Harmonic Content on R^ 

Considering sound pressure wave distortion due only to the second 
harmonic of p x (i.e., due to p 2 )» the nonlinear in-phase response factor 
was computed. The results are displayed graphically in figure 9 for 
various ratios of p 2 to p x . For this Mach number (M - 0.005), the 
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1 C * 4 10*3 10*2 10* 1 10 ° 

P' 

rms 

Figure 9. Response Factor Versus Root-Mean-Square 
Sound Pressure Level for Various Harmonic 
Contents; M = 0.005, ^ 6i = e 2 s 

T - T = 111.0 K° 

W 00 
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maximum response seems to occur for a nondimensionalized root-mean-square 
sound pressure level of 0.00631 and for P?/Pi approximately equal to 0.6. 
Figure 9 is similar in form to that obtained by Heidmann [3] for the 
combustion process; however, the peak response occurs at a lower value of 
^rms ^ or hot-wire than for the combustion process. This figure 
very clearly shows the effect of the distortion component p 2 on the 
response factor. 

Figure 10 was generated for the same conditions as figure 9 with 
the exception of the wire Mach number; figure 10 is for a zero through- 
flow velocity. This condition is of interest because it is the easiest 
to deal with experimentally and, as will be seen in figure 11, produces 
the maximum response. 

Influence of Mach Number on R » 

nc 

For a fixed pressure waveform (p 2 /Pi = 0.6) and for the value of 
P rms at which the peak response occurs, values of R^ were computed for 
various wire Mach numbers. The results are shown in figure 11. Figure 
12 shows the variation in with changes in for various wire Mach 
numbers. Generally speaking the lower the Mach number, the greater the 
response. It should also be noted that the value of the sound level 
at which peak response occurs increases with increasing Mach number. 

These results agree with those of Heidmann [3]. 

Influence of Harmonic Phase Angle on R^ 

The influence of *2 on the nonlinear response factor for the 
conditions of the peak R^ of figure 9, except for $ 2 » is shown in figure 
13. The response factor has its maximum at <f> 2 = 0 degrees. Figure 14 




Figure 10. Response Factor Versus Root-Mean-Square 
Sound Pressure Level for Various Harmonic 
Contents; M = 0, 4*1 = <j > 2 = = e 2 = 0, 








Response Factor Versus Phase Angle Between 

Fundamental Sound Pressure Component and 

Its Second Harmonic; M = 0.005, ^ = 0! = 

e 2 = 0, p 2 /Pi = 0.6, P’ = 0.00631 , 

T, - T = 111.0 K° rms 
w <» 



Figure 14. Response Factor Versus Root-Mean-Square 
Sound Pressure Level for Various Phase 
Angles Between Fundamental Sound Pressure 
Component and Its Harmonic; M = 0.005, 

(jjj = 6i = e 2 = 0, p 2 /Pi = 0.6, 

T - T = 111.0 K° 

W 00 
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shows the variation in R^ with root-mean-square sound pressure for 
various phase angles $2 between the fundamental and harmonic (distortion) 
component. Figure 15 is for the same conditions as figure 14 with the 
exception of the wire Mach number, which is zero in figure 15. Figures 
13, 14 and 15 show the very strong influence of the phase angle on the 
response. 


CONCLUSIONS 

The results obtained for the analog process — the output of a 
constant temperature hot-wire anemometer -- are qualitatively the same 
as those Heidmann obtained for vaporization limited combustion. There 
are small differences in the actual numbers involved, but this is to 
be expected since the equations for the two processes are slightly 
different. A scaling factor would be necessary to make a direct numerical 
comparison between the two, but this is unnecessary for the qualitative 
comparison done here. 

The addition of distortion to sinusoidal perturbations of the 
environment of the hot-wire can increase the response significantly. 

Using the stability criterion of Heidmann [3] -- response factors greater 
than some value between 0.8 and 1.0 denote unstable operation (the 
oscillations will grow with time) — the figures show that certain 
combinations of harmonic amplitude and phase angle will be unstable. 

The pure sinusoidal case, p 2 /Pi = 0, is always predicted to be stable. 
These results are examined experimentally in the next chapter. 




p' 

rms 

Figure 15. Response Factor Versus Root-Mean-Square 
Sound Pressure Level for Various Phase 
Angles Between Fundamental Sound Pressure 
Component and Its Harmonic; M = 0, 

<f>i = = = P 2 /P 1 = 0 . 6 , 

T - T = 111.0 K° 

w «> 



Chapter 3 


EXPERIMENTAL ANALYSIS 
INTRODUCTION 

The object of the experimental portion of the research is to 
determine the open-loop response of a constant-temperature hot-wire 
anemometer to sinusoidal and distorted sinusoidal acoustic oscillations, 
in the same way as was done analytically; and, through this, to gain 
insight into the closed-loop performance of the system. An experimental 
system was set up which would produce the same acoustic environment for 
the hot wire as was used in the analytical portion of the work (to the 
extent possible within the physical limitations of the system). With it 
the measurements required to determine the response could be made. The 
system consisted of four basic parts -- the test section itself, the 
means of generating the necessary acoustic oscillations, the means of 
making the necessary measurements to determine the response and the various 
types of equipment needed to calibrate and monitor the performance of all 
other parts of the system. 

Before arriving at the system which is described later in 
this chapter, the system went through many design modifications. Each 
of the four basic parts caused its own special problems. Some of the 
tests done on the interim system configurations resulted in important 
decisions concerning the scope of the research and the acoustic fields 
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used. These decisions had significant effects on the design of the 
equipment so they will be discussed first. 

No Through Flow 

It was originally planned that data would be taken with no through 
flow through the chamber (M = 0) and with various amounts of constant 
axial flow through the chamber, as was done analytically. A cylindrical 
test section was built with perforated end plates. To generate the flow, 
the test section was connected via a sheet metal cone and flexible 
tubing to a fan which drew air through the system. This system is shown 
in figure 16. 

When the chamber was tested to determine its resonant frequencies, 
it was found that the experimentally determined frequencies did not agree 
with those obtained from equation (A. 39). Further investigation showed 
that there was spatial distortion of the acoustic field due to the holes 
in the end plates. When the same resonance test was done with solid end 
plates on the chamber, the experimental and theoretical frequencies 
were virtually identical. The results of the tests are listed in 
table 1. Not all the predicted frequencies appeared due to the location 
of the microphone used to detect resonance. The microphone was at a 
pressure node for some of the resonant frequencies. Higher resonant 
frequencies also existed, but those shown are sufficient to establish 
a basis of comparison. 

Because of the distortion caused by the holes in the end plates, 
the analytical and experimental analyses would not be the same. Examina- 
tion of figure 12 showed that the no-flow condition is the one that is most 
sensitive to acoustic perturbations -- the value of the sound level at 
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Figure 16. Photograph of the First Experimental System 
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Table 1. Resonant Frequencies (77°F) 


Mode 

mnn z 

Theoretical 

(Hz) 

Holes 

(Hz) 

Experimental 

No Holes 
(Hz) 

001 

556 

583 

556 

100 

990 

1010 

990 

002 

1112 

1128 

1113 

101 

1135 

1192 

1135 

102 

1489 

1547 

1490 

200 

1642 

1661 

1651 

003 

1669 

1672 

1667 

201 

1734 

1765 

1750 

103 

1941 

— 

— 

202 

1983 

2020 

1981 

010 

2060 

2068 

2059 

on 

2134 

2151 

2134 

004 

2225 

2263 

2221 

300 

2259 

2284 

2257 

301 

2326 

2345 

2330 

012 

2341 

2374 

2339 

203 

2341 

2374 

2339 

104 

2435 

2491 

2432 

302 

2518 

2550 

2515 
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which the peak response occurs decreases with decreasing Mach number 
and the lower the Mach number, the greater the response factor. Hence, 
it was decided to restrict the experimental work to a truly closed ended 
cylinder with M = 0. It was felt that sufficient information could be 
obtained from just the no-flow situation to determine the validity of the 
analysis. 

Combination Mode 

Another major experimental difficulty, which was discussed briefly 
in Chapter 2, was caused by Heidmann's characterization of the distorted 
acoustic field as being the sum of the first spinning tangential mode and 
the second spinning tangential mode with 

f 2 DO = 2f 1 OD 

This is not true experimentally, as can be seen in table 1, or analytically, 
as indicated in equation (A. 38). Examination of this equation shows 
that for the first two pure tangential modes 


^200 a 20 lr 
flOO a 10 7r 


1.6604 


(3.1) 


and, therefore, f 2 oo is n °t equal to twice f 10 o* The ra tio of the two is 
always the same and it is determined by the first zeros of the Bessel 
functions 


and 


Ji ( a i o 7r ) ~ 0 

J2(cc 20 tt) = 0 


mn 


(3.2a) 

(3.2b) 


Values of a 


are given in Appendix B. 
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The ratio of the resonant frequency of a combination mode to that 
of the first tangential mode, however, is not only dependent on the zeros 
of a Bessel function but it is also dependent on the length and radius of 
the chamber, 

f " n z . (% n /R ) 2 * (yir 1 °' 5 

fioo («10/ R ) Z _ (3.3) 

The ratio of the two can be altered by changing the length or radius or 
both. Changing the radius changes f 100 and f mnn ; changing the length 
only affects f mnn . 

Examination of the acoustic equations for a (202) mode (equations 
(A. 48) with m = 2, n = 0, n z = 2) showed that although this mode has a 
longitudinal component of velocity, a longitudinal velocity node exists 
at the center of the chamber and that the equations for the pressure and 
tangential and radial velocities at the center have the same spatial 
relationships as the pure second tangential mode. Hence, the length 
of the chamber was adjusted so that 

2 02 ~ 2fioo (3.4) 

and the wire was positioned in the central cross section of the chamber. 
In this way the distorted acoustic oscillations used by Heidmann would be 
simulated as was shown in Chapter 2. 

APPARATUS 

A schematic diagram of the principal parts of the experimental 
system is shown in figure 17. Each of these parts is described, as 
is the additional equipment needed to calibrate and monitor the system. 
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Test Section 

The test section was originally designed to be part of a system 
which had a through-flow velocity and which would be excited by (100) 
and (200) pure spinning tangential modes of resonance. As was just 
discussed, this did not turn out to be the actual experimental condi- 
tions used. Rather than making a new test section, the original one was 
modified to meet the changing demands. Also, when any design modifica- 
tions were made, an effort was made to anticipate as many future needs 
as possible. Because of these things, the chamber has some features 
which are not needed for the present experiments, but which will be 
described for the sake of completness. 

Figure 18 is a photograph of the test section; it is an aluminum 
cylinder with a 20.49 cm (8.068 in.) I.D., 25.319 cm (9.968 in.) O.D. 
and length of 31.198 cm (12.283 in.). The cylinder is composed of two 
sections, one 12.70 cm (5.000 in.) long and the other 18.500 cm (7.283 
in.) long, which are mated by means of a lap joint containing an 0-ring 
seal. The longer section has provisions for mounting the hot-wire 
probe and the microphone from which the response factor is obtained, as 
well as other devices to examine the acoustic field. The shorter section 
has three threaded holes, for the mounting of acoustic drivers, located 
at its central transverse section, 90 and 180 degrees apart. When not 
in use, these holes are closed with threaded plugs which maintain the 
proper curved contour of the inside surface of the chamber. The 
cylinder was made in two sections so the angular position of the probe 
could be changed relative to the angular position of the drivers on the 
other section. 

The ends of the cylinder are closed with plexiglass end plates. 




Figure 18. Photograph of the Test Section 


<_n 



At one end are two 1.27 cm (0.50 in.) thick end plates. (The second one 
is for added structural rigidity.) Four 0.254-cm (0.10 in.) diameter 
holes are located at 0, 90, 180 and 270 degrees at r = 9.20 cm (3.62 
in.). These holes were made to permit investigation of the longitudinal 
variations of the pressure in the chamber. They are large enough to 
permit a 2.0 mm (0.079 in.) diameter microphone probe tube extension to 
pass through them. 

At the other end of the cylinder is a 5.08 cm (2.00 in.) thick 
end plate in which provision has been made for mounting three acoustic 
drivers and three microphones to monitor the sound field. The end 
plate must be thick to be able to support the heavy drivers without 
distorting the shape of the end plate. Plugs are available to close any 
holes not in use. Figure 19 is a photograph of the driver end plate with 
two drivers and three microphones mounted on it. The third driver hole 
is closed with a plug. 

Use of the (100) and (202) modes of vibration with the require- 
ment that 

^202 = 2fjoo (3.4) 

means the permissible values of length and diameter are functionally 
dependent, as is expressed in equation (3.3). To comply with this 
equation, the existing chamber either had to have a smaller radius, a 
greater length or modifications of both dimensions. To avoid having 
to make an entirely new chamber, the existing chamber was given a longer 
effective length by machining a disc 0.051 cm (0.020 in.) thick and of 
the chamber diameter out of each of the end plates. This resulted in 
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Figure 19. Photograph of an End Plate of the Test Section 
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the chamber having an effective overall length of 31.30 cm (12.323 in.) 
and meeting the requirement of equation (3.4). 

The cylinder and end plates have 0-ring seals between them and 
the whole system is clamped together with four long bolts. The chamber 
rests in a wooden cradle which has supports at both ends and at the 
lap joint. 

On the longer cylindrical section there are two probe access 
holes and over each of these is provision for the mounting of a probe 
holding and positioning device. One of the holes is located at the 
central cross-section of the cylinder and is 1.77 cm (0.50 in.) in 
diameter. This size permits the use of a 1.77 cm (0.50 in.) diameter 
microphone at this location, or smaller probes by means of appropriate 
sized plexiglass inserts. 

The second access hole is designed for use with a Disa right-angle 
miniature hot-wire probe holder. The hole itself is 1.905 cm (0.750 
in.) in diameter at the bottom and 3.81 cm (1.50 in.) in diameter at the 
top (refer to figure 17). The large size is necessary to be able to 
insert the right-angle probe holder. After insertion a plexiglass sleeve, 
made in two halves, is placed around the probe holder leaving only a 
0.474 cm (0.187 in.) diameter hole in which the 0.318 cm (0.125 in.) 
diameter probe holder can move. The center of the hole is located 7.00 cm 
(2.756 in.) from the central cross-section of the cylinder so that when 
the right-angle probe holder is being used the hot wire itself is 
located at the center of the cylinder and is on the same radial line as 
that which extends to the first probe site. Figure 20 is a photograph 
of the Disa right-angle probe holder, a hot-wire probe and the plexiglass 
insert. With suitable plexiglass inserts other probes can also be used 




Figure 20. Photograph of the Disa Right-Angle Miniature Probe Holder, A Hot-Wire 

Probe and the Plexiglass Insert m 

CO 
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at this location. The plexiglass inserts around the probes also act as 
electrical insulators between the probes and the test section and aid 
in preventing ground loops. 

Two probe holding positioners are available for use on the test 
section -- one manually adjusted and one electrically activated device. 
The manually adjusted holder has a vernier scale with which a probe can 
be positioned to within 0.25 mm (0.01 in.). It is shown in figure 21 
along with the inserts used to adapt it to the various probes. The 
electrically actuated holder (L. C. Smith BBR 24 180) has a remote 
control and readout device (L. C. Smith DI 3-R) associated with it — 
the probe location can be adjusted and the resulting position read from 
a digital voltmeter, with an accuracy of 0.65 mm (0.025 in.), at a 
location remote from the device itself. The photograph of figure 18 
shows this probe holder mounted on the test section. 

Acoustic Field Generation System 

The investigation requires that the combination of a (100) and 
(202) spinning acoustic field be generated and that the relative 
amplitudes and phases of these fields be adjustable with respect to one 
another. Since the frequency of resonance changes with temperature, and 
since chamber resonance studies are to be done, the frequency of the 
acoustic generation system also has to be adjustable. The number and 
location of the acoustic drivers is discussed first, then the means 
of generating the necessary oscillations to activate the drivers is 
explained. 

Driver location . As was discussed in Chapter 1, the acoustic 
drivers need to be located near pressure antinodes of the particular 
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Figure 21. Photograph of the Manually Adjusted Probe 
Holding and Positioning Device 



66 


mode desired. For the (100) and (200) modes, the choice includes any 
location for which r - R and 0 < z < L. For the (202) mode* the pressure 
antinodes are at r = R, but the longitudinal positions of the antinodes 
are restricted to 0, L/2 and L. The original location of the drivers — 
at the central cross-section of the shorter cylindrical section — was 
very close to a (202) pressure nodal position rather than an antinodal 
position. Thus the location of the drivers had to be changed. The center, 
z = L/2, was difficult to utilize since the microphone, the hot-wire probe 
holder and the cradle supporting the test section obstructed access. 
Elimination of the flow through the chamber (solid end plates) meant 
the ends of the chamber could be used as probe and driver sites; hence, 
the drivers were mounted on one of the end plates. 

In Appendix A, a spinning mode is shown to be the sum of two 
standing waves of equal amplitude but 90/m degrees out of phase with 
each other in angular position and 90 degrees out of phase in time. 

This indicates that to generate a spinning (100) mode, two drivers are 
needed separated by 90 degrees. To generate a spinning (202) mode, 
two drivers, separated by 45 degrees, are needed. Trigonometric manipula- 
tion of the equations shows that the same effect is achieved by putting 
the "second" driver at any pressure node of the acoustic field of the 
"first" driver. For the (202) mode this would be at 45, 135, 225 or 
315 degrees with respect to the first driver. 

Referring to figure 19, three driver sites are available — 1 
and 3, which are 90 degrees apart and therefore are suitable to drive 
the (100) mode, and 1 and 2, which are 135 degrees apart and are suitable 
to drive the (202) mode. In practice it was found that the (100) mode 
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could also be driven using 1 and 2; hence, only two drivers were actually 
used. The drivers are University Sound ID-75 driver units, rated at 75 
watts. Experience has shown that they cannot be used at this level when 
being driven by a pure tone, or a pure tone with one of its harmonics, 
without overheating. Most of the experiments were done with less than 
10 watts to each driver. 

Driver input signal generation. The system used to generate the 
signals required to activate the drivers was shown in figure 17. Figure 
22 shows the same system, only in more detail. A General Radio (GR) 
1309-A oscillator with a sinusoidal output is the basis of the sound 
generation system. The output frequency is adjustable from 10 Hz to 
100 kHz. The GR is tuned to the fundamental (100) frequency and its out- 
put is fed to the Wavetek 115 function generator, which uses it as a 
basis for generating a sine wave at the second harmonic (202) frequency. 
The Wavetek requires a signal of about 5 volts to be able to accurately 
frequency and phase lock to the GR. This is close to the maximum output 
of the GR. Because of this the amplitude adjustment of the GR cannot 
be used. Instead the amplitude is adjusted by means of an auxilliary 
voltage dividing potentiometer. When the frequency of the GR is changed, 
the frequency of the Wavetek also changes and maintains the proper 2 to 1 
ratio. 

The phase of the Wavetek signal can be adjusted relative to the 
GR signal. This phase adjustment is limited to about 80 degrees. The 
Wavetek has, in addition to its sinusoidal output, a negative sinusoidal 
output. By using this, an additional 80 degrees variation in phase can 
be obtained. When not in the frequency and phase lock mode, the output 





Figure 22 . Schematic Diagram of the System used to Generate the Acoustic Fields 
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frequency of the Wavetek is determined by the settings on its frequency 
control knobs, or is determined by the level of a D.C. input voltage. 

The amplitude adjustment of the Wavetek is not very sensitive. To 
obtain better control, a voltage dividing potentiometer was added. 

To activate driver 1, further fine adjustments of the amplitudes 
of each of the signals is made by means of two potentiometers; the two 
signals are then summed, fed to an Altec 1590-B power amplifier and then 
to driver 1. The summer is made from a Union Carbide UC4201 semiconductor 
integrated circuit operational amplifier, powered by a Kepco SC-36-1M 
regulated D.C. power supply. 

The signals to driver 2 each have to be phase shifted relative to 
the signals to driver 1. The amount of phase shift is different for 
each signal; hence, two phase shifters are needed. The fundamental 
frequency is shifted by means of a Lyrec TF 60 time-delay tape recorder 
operating in the loop mode; the harmonic frequency is shifted by means 
of an Acton 709 A broadband secondary phase standard. After the phase 
shifts, fine adjustments are made to the amplitudes of each of the signals 
by means of two potentiometers. The signals are summed, fed to a second 
Altec 1590 B power amplifier and then to driver 2. The second summer is 
identical to the first. When less complex acoustic fields are needed, 
parts of the system are simply not used. 

The photograph of figure 23 includes the GR and Wavetek function 
generators, potentiometer amplitude adjustors, the Acton phase shifter, 
and the Altec power supplies. The photograph of figure 24 includes the 
Lyrec tape recorder, which is used as a phase shifter, the amplifying and 
summing devices and the power supply for the summers. 
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Figure 23. Photograph which shows the B&K Frequency 
Analyzer, Fluke VAW Meter, Altec Power 
Amplifiers, GR Oscillator, Wavetek Function 
Generator, Voltage Dividing Potentiometers, 
Acton Phase Meter, Acton Phase Shifter and 
Stabiline Voltage Regulator. 
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Figure 24. Photograph which shows the Lyrec Time Delay 
Tape Recorder, B&K Level Recorder, Amplifying 
and Summing Devices and Associated Kepco 
D.C. Power Supply 
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The Lyrec recorder has a 300-inch tape loop adapter mounted on it. 
When operated in a loop, a blip appears in the output signal every time 
the splice in the tape which forms the loop passes over the read heads. 

To reduce the frequency of occurrence of the blip, the permissible 
length of the tape loop is increased by use of the adapter. Without it 
the recorder can handle only an 80-inch loop. 


Data Acquisition System 

The main property to be determined is the in-phase response factor 
It is a measure of the interaction between the anemometer output 
and the perturbations in the pressure and was defined as 
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No equipment is available which will determine this quantity directly; 
however, it is obtainable from the correlation coefficient, which can 
be determined experimentally. 

The correlation coefficient is defined as 

R = p- - - p! (3.5) 

rms rms 

The response factor is obtained from R as follows: 

From equation (2.30) 
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Equation 3.6 may be rewritten as 



Therefore , 


(3.7) 



Hence, to determine R^, an instrument which can convert the anemometer 

fluctuations E' and pressure fluctuations P' into the correlation 

coefficient is needed, as well as equipment to measure E' , P' „, the 

atmospheric pressure P and the time mean of the anemometer voltage F. 

The system used is shown in figure 17. Figure 25 shows it in more detail. 

The hot-wire probe used is a Disa 55F31 platinum-plated tungsten 

miniature hot-wire probe mounted on a Disa right-angle miniature probe 

support. Figure 20 shows both of these. The hot-wire is controlled 

by a Disa 55A01 constant temperature anemometer. The anemometer has 

two outputs -- a bridge output, E = E + E', and a turbulence output E'. 

The anemometer is shown in figure 26. 

The fluctuations in pressure P’ are sensed by a Briiel & Kjaer (B&K) 

4136 0.635 cm (0.25 in.) condenser microphone cartridge and B&K cathode 

follower. The power supply for this instrument is a B&K 2107 frequency 

analyzer, which can be seen in figure 23. The frequency analyzer is also 

able to measure the sound pressure level, from which P' , can be obtained, 

rms 

and perform frequency analyses of the microphone signal. 

The sound pressure level Lp is defined as. 
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Figure 25. Schematic Diagram of the System used to Obtain the Response Factor, R 
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Figure 26. Photograph which shows the Disa Random Signal 
Indicator and Correlator, Disa Constant 
Temperature Anemometer, TSI rms Voltmeter, 
Techni Lab Timer, H-P Integrating Voltmeter, 
H-P Oscilloscope, California Instruments 
Multimeter, L.C. Smith Probe Actuator, MKS' 
Baratron Pressure Transducer and Switch 



where 


Lp . 20 log 10 (P^ s /P ref ) 
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(3.9) 


P ref = 20 vN/m 2 


Solving equation (3.9) for P^, 


p • = p 

rms ref 


10 


Lp/20 


(3.10) 


(3.11) 


The outputs from the sound level meter and anemometer are fed to 

a Disa 55A06 random signal indicator and correlator, which determines 

the correlation coefficient R. The correlator is pictured in figure 26. 

A Thermo Systems Inc. (TSI) 1060 rms voltmeter is used to obtain 

E rms from E '' A Princo Fortin type barometer is used to obtain P. To 

determine E a Hewlett-Packard (HP) 2401C integrating digital voltmeter 

is used in conjunction with a Techni Lab Instruments II timer. They are 

interconnected so that when the timer is started, the voltmeter starts 

integrating with respect to time, E dt. When the timer is stopped, 

' o 

the voltmeter also stops. Dividing the integrated voltage by the time 
elapsed determines eT, i.e.. 


E = 


/. 


T 

E dt 


T 


(3.12) 


The TSI rms voltmeter, the HP integrating voltmeter and the timer are 
shown in figure 26. 

Calibration and Other Auxiliary Equipment 


To make a comparison between the analytical and experimental 
results, the equation relating the velocity to the voltage output from 
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the anemometer, equation (2.17), must be the equation which fits the 
particular hot wire being used in the experiment. To obtain this equation 
the hot wire must be calibrated. Also, the microphones used had to be 
calibrated and the acoustic fields generated had to be verified. In 
addition to the above, a resonant frequency analysis of the chamber was 
performed (as was mentioned at the beginning of this chapter); and the 
frequency of the various signals, the power to the drivers and the 
ambient temperature were monitored. Frequency spectrum analyses were 
done on both the anemometer and microphone outputs. Also, due to the 
intensity of the sound during operation, the test section had to be 
enclosed by an acoustic shield. 

Wire calibration equipment . The equipment used for the hot-wire 
calibration is shown photographical ly in figure 27 and schematically in 
figure 28. It consists of a Disa 55D41/42 low turbulence calibrating 
wind tunnel, an MKS Baratron 170M-5A pressure meter, a Neff 122 D.C. 
amplifier, a Meriam 34FB2 micromanometer, the anemometer, the integrating 
voltmeter-timer system described previously and a switch which allows 
either the output of the anemometer, or the output of the Neff amplifier, 
to be the input to the integrating voltmeter. 

The wind tunnel has provision for the mounting of a wire in the 
flow. The amount of flow is adjusted by a variable transformer which 
controls the tunnel blower motor. The velocity of flow is determined 
by measuring the static pressure difference between the inlet to the 
nozzle and the axial location of the hot wire. This pressure difference 
is sensed by the MKS Baratron pressure transducer, a capacitance mano- 
meter which converts the pressure difference to a voltage output. The 
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Figure 27. Photograph of the Equipment used for Hot-Wire Calibration 
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Figure 28. Schematic Diagram of the System used for Hot-Wire Calibration 
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voltage outputs for the velocities being considered (less than 3 m/sec) 
are small, so the Neff amplifier is used to increase them. The Baratron- 
Neff amplifier combination is calibrated using the Meriam micromanometer. 
The micromanometer is not sensitive enough to be able to accurately 
determine the pressure differences associated with velocities below 
about 1 m/sec. Hence, the micromanometer is used to calibrate the 
Baratron-Neff combination at higher velocities and the calibration, 
which is linear, is extrapolated to the lower velocities. 

Microphone calibration . To calibrate the microphones, a B&K 
4220 pistonphone is used. The pistonphone produces a sound pressure 
level of 124.0 dB (+ 0.2 dB } plus a correction factor determined by the 
barometric pressure. A barometer is supplied with the pistonphone 
which converts the barometric pressure reading to the necessary correction 
factor. Adaptors are available to fit the opening of the pistonphone 
to various microphones. Figure 29 shows the pistonphone being used to 
calibrate a microphone. The barometer and adaptors are also shown. 

Acoustic field verification . Prior to taking data, the presence of 
the desired acoustic fields had to be verified. Initially, rather 
extensive investigations of the acoustic fields were performed. Once 
the ability of the equipment to generate the necessary acoustic fields 
was demonstrated, the procedure for setting up the acoustic fields was 
followed and the resulting sound field was used for actual data taking 
without further verification. 

From equation (A. 48a) the variations in the pressure in the 
chamber with respect to both space and time can be determined. To 
verify the existence of a particular mode of vibration, these analtyical 
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Figure 29. Photograph of the Microphone Calibration 
Equipment 
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predictions of the pressure profiles, and phase relationships of the 
pressure in various parts of the chamber, are compared to the experi- 
mentally determined ones. 

To measure a radial pressure profile a B&K UA0040 probe micro- 
phone kit was used to make a 2 mm (0.079 in.) O.D., 22.86 cm (9.0 in.) 
long microphone probe which was attached to the B&K 4136/2615 micro- 
phone (microphone 3). A plexiglass adaptor was made for the central 
probe hole of the chamber to reduce the size of the hole to 2.36 mm 
(0.093 in.). The microphone probe was clamped into the manually 
adjusted probe holder (figure 21), which was mounted over the probe 
hole and used to position the tip of the probe at various radial 
locations in the chamber. The amplitude and the phase of the pressure 
was measured at each location. The phase measurement was accomplished 
by comparing the appropriate oscillator signal to the microphone signal, 
using a Acton 320- AB phase meter. The same microphone- probe adaptor 
combination was used in a similar fashion to investigate the longitudinal 
variation in the amplitude and phase of the pressure. 

To examine the tangential variation in amplitude and phase, three 
additional microphones were used which were mounted in the same end 
plate as the acoustic drivers. These microphones are numbered as shown 
in figure 19. Microphone 1 is a B&K 4133 1.27 cm (0.50 in.) condenser 
microphone cartridge with a B&K 2615 cathode follower. It is powered 
by a B&K 2603 microphone amplifier. Microphone 2 is a B&K 4134 1.27 cm 
(0.50 in.) condenser microphone cartridge with a B&K 2615 cathode 
follower. It is powered by a B&K 2801 microphone power supply. Micro- 
phone 4 is a B&K 4135 0.635 cm (0.25 in.) condenser microphone cartridge 
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with a B&K 2615 cathode follower. It is powered by a B&K 2803 two- 
channel power supply. To determine the sound pressure level for any of 
these, they are connected to the B&K 2107 frequency analyzer. Each 
microphone is calibrated using the pistonphone. The three microphone 
power supplies can be seen in the photograph of figure 18. 

When both the (100) and (202) acoustic fields exist in the chamber 
at the same time, each may be analyzed separately by using the frequency 
analysis mode of the B&K 2107 frequency analyzer. In this mode the 
particular frequency of interest is chosen by the setting of a frequency 
range dial and frequency tuning scale. Six different bandwidths are 
available. 

Resonant frequency analysis . When the air in a chamber is excited 
by an acoustic driver, the resulting sound pressure level, for a fixed 
power input, depends on the frequency of the driver input. The sound 
pressure level peaks at the resonant frequencies of the chamber. This 
fact was used to perform the resonant frequency analysis mentioned at 
the beginning of the chapter. 

The system shown in figure 30 was used to determine the resonant 
frequencies of the chamber. The acoustic driver is activated by the 
Wavetek function generator and the Altec power amplifier. The output 
frequency of the Wavetek is determined by the input voltage level to 
it. This voltage is supplies by a Hewlett-Packard 3310A function 
generator. The HP is set to deliver a very low frequency ramp function. 
The frequency used was about 2 x 10 -4 Hz, which means one cycle would 
take almost an hour and a half to complete. In practice the total 
cycle was not used; the system was stopped when it was felt a sufficient 
range of frequencies had been examined. The amplitude was adjusted so 




Figure 30. Schematic Diagram of the System used to Determine the Resonant 
Frequencies of the Chamber 


00 

-p* 




85 


that frequencies from about 100 Hz to greater than 4000 Hz were available 
The frequency being produced by the Wavetek was measured by a California 
Instruments series 8300 digital multimeter. The HP function generator 
and multimeter can be seen in figure 26. 

The microphone used to detect the sound pressure level is powered 
by the B&K 2107 frequency analyzer. The output from the 2107 is fed 
to a B&K 2305 level recorder, which plots the sound pressure level as 
a function of time. Since the frequency is dependent on time, the 
abscissa is also a measure of the frequency. To calibrate the abscissa, 
a mark was made on the chart at every 100 Hz interval, as indicated by 
the multimeter. The frequency varied linearly with time. The resonant 
frequencies are the frequencies at which the peaks occur on the chart. 

The level recorder is shown in figure 24. 

Miscellaneous equipment . The power to the drivers, as well as 
the voltage and current, was monitored by a John Fluke 102 VAW meter. 

A switch was incorporated into the circuits of the drivers which enabled 
the VAW meter to be switched to either driver. The VAW meter is 
included in the photograph of figure 23. 

An additional frequency meter was also used, a General Radio 
1191-B counter. With it frequencies can be measured to within 0.1 Hz. 

It can be seen in figure 26. 

To better understand the nature of the anemometer and microphone 
signals, spectrum analyses were performed on them. To do this an EMR- 
Schlumberger 1510 spectrum analyzer was used. A Houston Instruments 
Omni graphic 2000 X-Y recorder was interconnected with the spectrum 
analyzer and was used to record the spectrum obtained. The analyzer 
and plotter are shown in figure 31. 




Figure 31. Photograph which shows the EMR-Schlumberger Spectrum Analyzer, Houston 
Instruments X-Y Recorder, Tektronix Oscilloscope and Polaroid 

Oscilloscope Camera g 
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To study the phase relationships between the pressure and anemo- 
meter signals, and the shape of these signals, a Tektronix 502-A dual- 
beam oscilloscope was used. A Polaroid oscilloscope camera (Tektronix 
C-27) was mounted on the oscilloscope and pictures were made of the 
signals. The phase difference between the signals could then be 
determined from measurements of the distance between the peaks in the 
pictures. The oscilloscope and camera can be seen in figure 31. An 
additional oscilloscope was also used to monitor the wave shapes of the 
various signals — a Hewlett-Packard 1200B dual -trace oscilloscope 
which can be seen in figure 26. 

A Stabiline EM 10364 voltage regulator was used to supply 
regulated A.C. power to all of the electronic equipment. 

When the system was in operation, it was found that the sound 
levels were too high for the personal safety or comfort of those 
operating the system, and were annoying to others working in the same 
building. Because of this, an acoustic shield was built which housed 
the test section itself. Most of the electronic equipment was kept on 
the outside of the shield. The shield is a structure 12-feet long, 
8-feet wide and 8-feet high. The walls and roof are formed by a 2 x 4 
framework covered with two layers of 3/8- inch sheetrock on the inside 
and two layers of 3/8- inch sheetrock on the outside. The roof is sealed 
to the top of the walls, and the bottom of the walls are sealed to the 
floor, with rubber weather stripping. The door is also sealed all 
around with rubber weather stripping. The inside is furnished with a 
fluorescent ceiling light and power receptacles. The cables connecting 
the chamber and the electronic equipment pass through a 2.5-inch 
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diameter hole in one of the walls. When in operation, the sounds being 
generated can still be heard outside the shield, but at a much reduced 
level. The acoustic shield is shown in figure 32. 

An overall view of the electronic equipment is shown in figure 33. 

PROCEDURES 

All the electronic equipment was kept on at all times. Period- 
ically, the various instruments were adjusted according to the 
manufacturers' instructions. 

Six procedures will be described -- microphone calibration, 
chamber resonant frequency analysis, acoustic field generation and 
verification, hot-wire calibration and response factor determination. 

Microphone Calibration 

To calibrate a microphone the following procedure is used: 

® The microphone is connected to its power supply and to the 
sound pressure level meter. 

© The microphone is inserted into the pistonphone's 
cavity (using an adaptor is necessary). 

© The pistonphone is turned on and the sound pressure level 
on the meter is determined. The pistonphone is turned off. 

® The sound pressure level in the cavity is calculated by 
measuring the barometric pressure, determining the cavity 
correction factor from the barometer and adding this 
factor to 124.0 dB. 

® The difference between the sound pressure level produced in 
the pistonphone cavity and the meter reading is determined. 
This difference is the microphone calibration factor; it 




Figure 32. The Acoustic Shield 
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Figure 33. Overall View of the Electronic Equipment 
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must be added to the meter reading to obtain the true sound 
pressure level. 

For the 2107 frequency analyzer the calibration factor depends on 
whether or not the instrument is being used in the frequency analysis 
mode. The microphones must be calibrated for both conditions. 

The measured phase of a signal varies slightly depending on which 
microphone is used to detect it. Since the phase difference between 
various locations in the chamber is an important part of verifying that 
the correct acoustic field exists in the chamber, the microphones must 
be calibrated for phase relative to one another. This is done as follows: 

© A standing resonant acoustic field is set up in the chamber. 

© The microphone which is chosen as the "standard" is placed 
at a pressure antinodal position. 

© The phase difference between the oscillator used to generate 
the acoustic field and the microphone is measured using the 
Acton phase meter. 

® The "standard" microphone is removed and the other micro- 
phones are each successively placed at the same location 
and their phases are also measured. 

® The difference between the "standard" phase measurement 
and each of the other microphones' phase measurement is then 
determined. The resulting numbers are the phase calibration 
factors for the microphones. This calibration must be 
done for each freqency used. 

Chamber Resonant Frequency Analysis 

To perform the chamber resonant frequency analysis, the system 
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is set up as shown in figure 29, with the exception that the HP and 
Wavetek function generators are not interconnected. The following 
procedure is then followed: 

® The HP function generator is set for a ramp function 
output at a frequency of about 2 x 10~ 4 Hz and an 
amplitude of about 5 volts. 

© The Wavetek function generator is scanned manually 
from about 500 Hz to about 3000 Hz to determine the 
maximum sound pressure level to be expected. 

® The B&K recorder is scaled so the maximum sound pressure 
level will not go off scale and an appropriate paper speed 
is set — about 0.3 mm/sec (0.71 in. /min). 

• When the HP function generator reaches the beginning of 
a cycle, it is connected to the Wavetek. This initiates 
the frequency scan. 

© The B&K recorder is started. 

© As each 100 Hz interval is reached, as indicated by the 
multimeter, the chart paper on the recorder is marked. 

© When a sufficient range of frequencies has been scanned, 
the recorder is turned off and the HP is disconnected. 

# The frequencies at which the peaks occur on the chart are 
measured. These are the resonant frequencies of the chamber. 

Acoustic Field Generation and Verification 

The desired acoustic field is a (100) spinning mode, distorted 
by a (202) spinning mode. These two modes are set up on the following 


basis: 
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@ A spinning wave is the sum of two standing waves. 

@ The sound pressure level of a spinning mode, in a transverse 


section, at a particular value of r, is constant for any 
value of e. 

9 For a (100) spinning mode the phase difference between the 
pressure at e = 0j and the pressure at e - e 2 is equal to 
(01 - 0 2 ). 

© For a (202) spinning mode, in a transverse section, the 
phase difference between the pressure at e = and the 
pressure at e = 0 2 is equal to 2(0! - e 2 ). 

To generate the acoustic field, the system shown schematically 
in figure 22 is used. The drivers and microphones are located as shown 
in figure 19. The microphones all have different power supplies, as 
was described previously, but they use the same B&K 2107 frequency 
analyzer for sound pressure level determinations. The frequency analysis 
mode is used to isolate the pressure associated with the (100) mode 
from that associated with the (202) mode. The microphones are calibrated 
for sound pressure level and phase. Prior to setting up the acoustic 
field, the inputs to the summers and the amplifiers are not connected. 

All other connections are as indicated in figure 22. The power amplifiers' 
controls are set at 8. The following is the procedure used to set up the 
acoustic field: 

1. The amplitude adjust of the GR is turned up to an 

"appropriate" level. (An appropriate level depends on 
the sound pressure level desired. A sound pressure level 
of 155 dB or more would require close to the maximum 
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output. A sound pressure level of 140 dB would only 
require about a third of that.) 

2. The GR output is connected to amplifier and summer 1 (AS1). 

The output of AS1 is connected to power amplifier 1 (PA1). 

This sets up a standing acoustic field in the chamber. 

3. Microphone 1 (Ml) is connected to the B&K 2107 and the 

frequency of the GR oscillator is tuned to the (100) 

resonant frequency. This frequency varies slightly 
with the temperature, but it is approximately 990 Hz. 

Tuning is accomplished by varying the frequency from 

990 Hz until the maximum sound pressure level is obtained. 

4. The B&K 2107 is turned to the frequency analysis mode and 

set to the (100) frequency. 

5. The sound pressure level is adjusted to the approximate level 
desired for the (100) mode by use of the appropriate potentio- 
meter of AST. 

6. The input to PA1 is disconnected. The phase shifted GR 
signal is connected to AS2 and the input to PA2 is connected. 

7. M2 is connected to the B&K 2107 and the sound pressure level 
is adjusted to the approximate (100) level desired, by use 
of the appropriate potentiometer of AS2. 

8. The input to PA1 is reconnected. Two standing waves of 
approximately the same amplitude now exist in the chamber. 
Next, the proper phase relationship between the two must 
be established. 

9. Ml is connected to the B&K 2107. The exact (100) sound 
pressure level desired is set by adjusting the appropriate 
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potentiometer of AST . The phase between the output of 
the GR and Ml is determined using the Acton phase meter. 

10. M2 is connected to the B&K 2107. The phase between the GR 
and M2 is measured using the phase meter. The Lyrec phase 
shifter is adjusted until the difference in the phases 
measured in steps 9 and 10 is 135 degrees. The exact (100) 
sound pressure level desired is set by adjusting the 
appropriate potentiometer of AS2. 

11. Steps 9 and 10 are repeated until the sound pressure levels 
of Ml and M2 are the same and the signals are 135 degrees out 
of phase. A (100) spinning acoustic mode should now exist in 
the chamber. Next, the (202) spinning mode is set up. 

12. The amplitude adjust of the Wavetek is turned up to an 
appropriate level. 

13. The Wavetek output is connected to AS! . This sets up a 
standing (202) acoustic field in the chamber. It is at 
resonance because the Wavetek frequency is twice the GR 
frequency and the (202) resonant frequency is twice the 
(100) resonant frequency. Tuning for the (100) mode tunes 
for the (202) mode. 

14. The B&K 2107 is set to the (202) frequency. (It is still 
in the frequency analysis mode.) 

15. The sound pressure level is adjusted to the approximate level 
desired for the (202) mode by use of the appropriate 
potentiometer of AS1 . 

The phase shifted Wavetek signal is connected to AS2. 


16 . 
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17. M2 is connected to the B&K 2107 and the sound pressure level 

is adjusted to the approximate level desired for the (202) 

mode by use of the appropriate potentiometer of AS2. 

18. Ml is connected to the B&K 2107. The exact (202) mode 
sound pressure level desired is set by adjusting the 
appropriate potentiometer of AS1. The phase between the 
output of the Wavetek and Ml is determined using the phase 
meter. 

19. M2 is connected to the B&K 2107. The phase between the 

Wavetek and M2 is measured. The Acton phase shifter is 

adjusted until the difference in the phases measured in 

steps 18 and 19 is 270 degrees. The exact (202) mode 
sound pressure level desired is set by adjusting the 
appropriate potentiometer of AS2. 

20. Steps 18 and 19 are repeated until the sound pressure levels 
of Ml and M2 are the same and they are 270 degrees out of 
phase. A (202) spinning acoustic mode should now exist 

in the chamber along with the (100) spinning acoustic 
mode. 

21. Minor adjustments in the levels of the two modes can be 
made using the GR and Wavetek amplitude adjustment 
potentiometers. 

22. The phase for the (202) mode with respect to the (100) mode 
can be changed by varying the phase adjustment knob of the 
Wavetek. 

Before proceeding to studies of the response factors using this 
method to set up the acoustic field, more extensive investigations were 
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performed to verify that the desired acoustic field had been set up. The 
variation in the level and phase of the sound pressure with respect to e, 
r and z is predicted by equation (A. 48a). Measurements were made of 
the sound field to see if it matched the theory. 

When setting up the sound field for verification purposes, two 
additional steps are added to the procedure: 

11a. M4 is connected to the B&K 2107 and the sound pressure 
level and phase between the Gr and M4 is measured. 

20a. M4 is connected to the B&K 2107 and the sound pressure 
level and phase between the Wavetek and M4 is measured. 

These two steps are an additional check on the variation of sound 
pressure level and phase with respect to e. 

The procedure for checking the variation in level and phase with 
respect to r consists of the following steps: 

1. Microphone 3 in combination with the probe tube extension is 
positioned over the central probe access hole using the 
manually adjusted probe holding and positioning device. 

M3 is connected to the B&K 2107 frequency analyzer. 

2. The sound field is set up in the chamber using the method 
described previously. 

3. With the B&K 2107 in the frequency analysis mode and set to 
the (100) frequency, values of the sound pressure level and 

the phase with respect to the GR function generator are measured 
at various radial locations in the chamber. 

4. With the B&K 2107 in the frequency analysis mode and set to 
the (202) frequency, values of the sound pressure level and 
the phase with respect to the Wavetek function generator 
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are measured at various radial locations in the chamber. 

The procedure for checking the variations in level and phase with 
respect to z is the same as that with respect to r, except that the 
microphone probe and holder are mounted on an end plate, 

Hot Wire Calibration 

The wire calibration is divided into two parts data taking 
and data analysis. 

Data taking . The following is the procedure used for taking the 
data from which to calibrate a hot wire: 

1. The calibration equipment is set up as shown in figure 28. 

The variable transformer is turned to the off position. 

The hot wire is mounted in the wind tunnel in the orienta- 
tion that will be used in the experimental system. 

2. The operating temperature of the hot wire is chosen 

and the wire operating resistance associated with this is 
set on the anemometer. The anemometer is turned to 
"operate. " 

3. The micromanometer is set to zero, 

4. An appropriate amount of amplification is set on the 
Neff amplifier. 

5. The barometric pressure and ambient temperature are 
measured. 

6. The switch is set to "E.” The timer and the HP integrating 
voltmeter are reset to zero. The timer is turned on, 
which starts the integration. After a minimum elapse 

time of 100 seconds, the timer is turned off, which stops 
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T 

the integration, and the resulting values of J E dt and T 
are recorded. 

7. The switch is set to "H. " The timer and the voltmeter are 
reset to zero. The timer is turned on, which starts the 
integration. After a minimum elapse time of 100 seconds 
the timer is turned off, which stops the integration, and 
the resulting values of /J H dt and T are recorded. 

8. The setting on the variable transformer is increased. 

This places a flow across the wire. Steps 6 and 7 are 
repeated. 

9. Step 8 is repeated until a sufficiently high flow rate 

has been established so that readings of the pressure drop 
across the nozzle can be read on the Meriam micromanometer; 
a micromanometer reading of aP is taken after steps 6 and 7 
are completed. Care must be taken that sufficient time 
elapses between increasing the flow rate and the taking of 
the data so that the system has had time to reach steady state. 

10. When a sufficiently high velocity has been reached, the system 
is shut down and the data are analyzed. 

Data analysis . The hot wire calibration equation is 

E = [A + 8 V N ] 0 - 5 (2.18) 

What follows is the procedure used to fit the data to an equation of this 
form; the constants A, 8 and N are to be determined: 

1. Each value of /J E dt and /J H dt is divided by the value 
of T associated with it. This determines time-average 
value of E and H for each flow rate. 
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2. The H data, which are in volts, must be converted to 
equivalent pressures. This is accomplished by performing 
a least-squares curve fit of the H versus AP data; the 
relationship is linear. When the equation is obtained, 
the values of H are converted to aP. 

3. The velocity V which corresponds to each value of aP 
is found using the following equation: 

V 2 = [^r] R a T[1 " (1 " Ap / p ) {Y_1/Y) ] (3.13) 

This equation is taken from the instruction manual for the 
calibrating wind tunnel [20]. 

4. Equation (2.18) can be put in the form 

In (E 2 - A) = In 8 + N In V (3.14) 

Letting 

In (E 2 - A) = Y (3.15a) 

In B = b (3.15b) 

and 

?n V = X (3.15c) 

equation (3.14) becomes 

Y = b + NX (3.16) 

A value is assumed for A. Using this assumed value, the data 
are converted to Y and X and a computer program is used to fit 
the resulting values to equation (2.18). The program 
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determines a value for b and N and also a value for the 
"goodness-of-fit" of the equation to the data. 

5. Other values of A are assumed and the procedure in step 
4 is repeated for each. 

6. A plot of the goodness-of-fit versus the assumed value 
of A is made. The value of A where the goodness-of-fit 
reaches a maximum is determined. This value of A, along 
with the corresponding values of b and N, determine the 
wire calibration equation constants A, 8 and N. 

By knowing the temperature difference at which the calibration 
was performed, together with various other wire and flow field 
characteristics, equation (2.17) can be obtained from equation (2.18). 
Equation (2.17) is then used for the analysis done in Chapter 2. A more 
detailed discussion of the hot-wire calibration procedure can be found 
in reference [20]. 

Response Factors 

To obtain the data necessary to calculate response factors, 
the experimental system shown schematically in figure 17 was used. 

Both the hot-wire and microphone probes should be at the same 
location. Since this is impossible, some compromise had to be made. 
Also, the interference of the probes with the sound field should be 
kept at a minimum. 

The rate of change of P^ ms with respect to radial position 
decreases as r increases and approaches zero as r approaches R. Because 
of this, and to prevent the microphone from affecting the sound field, 
the microphone used to obtain the response factors (B&K 4136, microphone 
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3) was mounted flush with the inside surface of the test section. (As 
was discussed previously, the longitudinal location of the probes, 
z = L/2, was determined by the characteristics of the (202) acoustic 
field. ) 

The hot wire was positioned directly under the microphone at 
r/R = 0.875 and was oriented vertically (coincident with the radius 
vector to the microphone). 

The quantity 4> 2 - 4>j is the phase difference between the funda- 
mental and distortion component of the acoustic oscillations using the 
terminology Heidmann developed. The relationship between this terminology 
and that which was used in the solution of the inviscid wave equation 
was explained in Chapter 2. The phase difference can be measured using 
any arbitrary reference point; therefore, <j>j will be assumed to be 
zero and, henceforth, the phase difference <t> 2 - <j>i will be equal to 
4 > 2 . 

Two procedures were used to obtain response factors. The first 
procedure generated the variation in with respect to for <f> 2 

equal to zero and 180 degrees and for a constant value of p 2 /Pi* The 
second procedure generated the variation in R„ with respect to <j> 2 
when P^, and p 2 /Pi are held constant. 

as a function of P^ ms - With the system set up as illustrated 

in figure 17, the procedure to determine R . as a function of P* is as 

nt rms 

follows: 

1. A value of p 2 /Pj is chosen. From this the corresponding 
difference between the sound pressure levels ( LPi 0 o ~ Lp 202 ) 


is determined. 
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2. All the microphones are calibrated. 

3. The atmospheric pressure and temperature are recorded. 

4. The anemometer is switched from its stand-by to its 
operating condition. 

5. The signal from the Wavetek is taken from the positive 
sinusoidal output. The Wavetek's phase shift knob is set 

at its minimum position (counter clockwise). At this setting 
it is least likely to become unlocked while manipulations 
are being made to set up the acoustic field. 

6. With the B&K 2107 frequency analyzer in the frequency 
analysis mode and using microphones 1 and 2, an acoustic 
field is set up which has approximately the required value 
of (LPioo “ Lp 202 ) and which has the maximum overall sound 
pressure level obtainable without overheating the drivers 
(the power to the drivers is limited to less than 10 watts). 

7. Using microphone 3 (with the B&K 2107 still in the frequency 
analysis mode) and the external amplitude adjustment 
potentiometers of the GR and Wavetek function generators, 
the sound pressure levels are adjusted to the exact value 

of (LPioo - LP 202 ) desired. 

8. With the B&K 2107 no longer in the frequency analysis mode 
(all frequencies now included in the microphone signal) the 
phase shift of the Wavetek is adjusted until the trace of 
the pressure signal, as seen on an oscilloscope, has "even 
peaks" ( f^X/^X ). This is the waveshape predicted when 
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<t> 2 is equal to zero degrees. 1 

9. The value of (Lp 100 - Lp 20 2 ) is rechecked. 

10. The overall sound pressure level is read from the meter 
on the B&K 2107. 

11. The timer-integrating voltmeter combination is started. 

12. The value of E' „ is read from the TSI rms voltmeter. 

rms 

13. The necessary adjustments are made to the correlator 
controls and the value of R is read from its meter and 
recorded. 

14. The timer-integrating voltmeter system is stopped and the 
resulting values of /J E dt and T are recorded. Sufficient 
information has now been obtained to determine a value of R^. 

15. The Tektronix 502A oscilloscope and camera are used to make 
a photographic record of the microphone and anemometer 
signals. 

16. The EMR-Shlumberger spectrum analyzer is used to obtain the 
frequency spectrum of the microphone and anemometer signals. 

A permanent record is made of these spectrum analyses using 
the Texas Instruments X-Y plotter. (It is not necessary to 
carry out steps 15 and 16 every time since they are not 
necessary to determine the value of R nV 


x No method was available to measure the absolute difference in 
phase <t> z between the (100) and (202) acoustic fields; only the change 
in phase could be measured. Hence, the even peaks condition associated 
with <t> 2 ~ 0 was used as a reference, and other phase angles were 
determined by adding or subtracting, whichever was appropriate, the 
change in <j> 2 to it. The method used to measure the change in <t> 2 
is described in the procedural description following this one. 
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17. The output from the Wavetek is changed from the positive 
to the negative sine and the phase adjustment knob is used 
to set up "even peaks" again. The peaks are now on the 
bottom instead of on the top ( VAJ ), which corresponds 
to <j >2 = 180 degrees. 

18. Step 7 is repeated. 

19. Steps 10 through 14, 15 or 16 are repeated. 

20. Steps 5 through 19 are repeated, but at a lower overall 
sound pressure level. 

21. Step 20 is continually repeated until the overall sound 
pressure level is so low that meaningful measurements can 
no longer be made. 

as a function of <p 2 - The procedure used to obtain R^ as a 
function of <f> 2 is as follows: 

1. With the system set up as illustrated in figure 17, steps 

1 through 7 of the previous procedure are followed with the 
exception that the overall sound pressure level is not 
necessarily set to the maximum obtainable. Usually it is 
more desirable to use the sound pressure level which 
corresponds to the maximum value of R^ as determined by 
previous data or analysis. 

2. The phase angle between the outputs of the GR and Wavetek 
function generators is measured using the Acton phase meter. 
This establishes a relative phase difference between the 
(100) and (202) acoustic fields. It is not an absolute 
phase difference since the phase changes that occur between 
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the function generators and the actual sound fields is not 
known. 

3. Steps 10 through 14, 15 or 16 of the previous 
procedure are followed. 

4. The relative phase between the GR and Wavetek function 
generators is changed by a few degrees using the phase 
adjustment knob of the Wavetek. The new relative phase 
difference is measured using the Acton phase meter. The 
change in phase is the difference between the previous phase 
measurement and this measurement. 

5. Steps 3 and 4 are repeated until the Wavetek is no longer 
able to remain locked to the GR. During this step the phase 
angle reading at which the even peaks condition occurs is 
recorded. 

6. The phase adjustment knob of the Wavetek is returned to its 
counterclockwise position, the output of the Wavetek is 
changed from the positive to the negative sine and the 
frequency and phase lock condition of the Wavetek is 
reestablished. 

7. Steps 2 through 5 are repeated. 

When these two procedures have been completed, the response 
factors are calculated from the data generated using equation (3.8) 



P 

E 


(3.8) 
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RESULTS AND DISCUSSION 

The results of a hot-wire calibration and the sound field 
verification procedures are presented and discussed. Data were taken 
to determine the microphones' and anemometer's operational characteristics 
certain of these characteristics are depicted. Finally, the response 
factors themselves are presented. 

The results of the chamber resonant frequency analysis have 
already been presented in the introduction of this chapter. 

Hot-Wire Calibration 

Two hot wires were calibrated during the course of the research. 

The resulting calibration equations are similar. Figure 34 contains the 
data used to determine a calibration equation and the equation itself. 

For this case the calibration constants of equation (2.18) are: 

A = 11.909 
8 = 8.0995 
N = 0.444 

with a goodness-of-fit of 0.99988. Combining these values with the 
various wire and flow field properties, the constants of the anemometer 
equation, equation (2.17), are found to be 

A = 0.525 
B = 0.6729 

The equations and data used to obtain A and B are contained in 
Appendix C. These constants were used in the analytical portion of 
the research. 



Anemometer Voltage, E (volts) 



Figure 34. Hot-Wire Calibration Data and Resulting Equation 
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The anemometer equations were developed assuming forced convection 
heat transfer from the hot wire to the surrounding cooler air. At a 
velocity of zero, or very close to zero, free convection heat transfer 
dominates and the equations are not valid. Because of this, the 
anemometer output at zero velocity is not a valid data point for obtaining 
the wire calibration constants and, therefore, it is not used in the 
least-squares curve fitting procedure. 

To account for the change in heat transfer mode in the vicinity 
of zero velocity, the calibration equation is considered to consist of 
two regions -- a very short free convection dominated region in which 
the anemometer output E is assumed to remain constant at its zero 
velocity value and a forced convection dominated region in which the 
equation for E is that determined by the calibration procedure described. 
The boundary of these two regions is the intersection of the constant E 
region and the variable E region. These two regions are shown in 
figure 34. 

Sound Field Verification 

Verification of the sound field depends on comparison of the 
analytically predicted sound pressure in the chamber with that obtained 
experimentally. The expected analytical distributions were discussed 
in Chapter 1 . 

Figure 35 shows the variation in the normalized values of P' 

ms 

across a diameter of the chamber, at the central cross section, for the 
(100) mode. Information regarding the phase of the pressure along the 
diameter with respect to the pressure at r = R is also shown in this 
figure. Figure 36 is similar to figure 35, but it is for the (202) mode. 




Phase 

' / f p ' ) 

rms 7 rms' max (degrees) 



r/R 

Figure 36. Sound Pressure and Phase Variation of the Second Tangential - Second Longitudinal 
Mode Along a Diameter at the Central Cross Section of the Chamber 
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Figure 37 shows the variation in the normalized value of P^ s and 
phase with respect to longitudinal position, at r/R - 0.90, for the (202) 
mode. The same data were taken for the (100) mode; however, no variation 
in either pressure of phase is predicted for the (100) mode and the 
data indicated only a slight decrease in P^ ms near the end plate. 

For the spinning situation, there should be no tangential 
variation of pressure amplitude for a constant value of r. The 
difference in the values of sound pressure level of the (100) mode, as 
measured by the three microphones mounted in one of the end plates, 
varied from 0.2 dB to 0.7 dB. The difference in the values of sound 
pressure level of the (202) mode, as measured by the three microphones 
in the end plate, varied from 0.3 dB to 1.3 dB. 

The above errors in the sound pressure level are related to the 
errors in the value of the phase angle measurements between the 
microphones. As was explained earlier, for a (100) spinning mode, 
the phase difference between the pressure at e = and the pressure at 
6 = e 2 is equal to (ej - e 2 ). The error associated with the phase 
measurements between the three microphones in the end plate for the 
(100) mode, ranged from 2 to 8 degrees. For a (202) spinning mode, the 
phase difference between the pressure at e = e r and the pressure at 
9 - 0 2 is equal to 2(e x - e 2 ). The error associated with the phase 
measurements between the three microphones in the end plate, for the 
(202) mode, ranged from 1 to 8 degrees. 

The larger errors in the sound pressure levels are associated 
with the smaller errors in the phase measurements, and vice versa. 

The errors are in part caused by the dimensions and locations of the 
microphones. Microphones 1 and 2 have a 1.27 cm (0.50 in.) diameter 
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Figure 37. Sound Pressure and Phase Variation of the Second Tangential - Second Longitudinal 
Mode with Respect to Longitudinal Position at r/R = 0.90. 



and are mounted so their circumference is tangent to the circumference 
of the cylinder. The third microphone, microphone 4, has a 0.635 cm 
(0.25 in.) diameter and is mounted so its circumference is 0.317 cm 
(0.125 in.) from the circumference of the cylinder. The measurements 
are being made over a finite area, not a point. In addition, the 
acoustic field is being driven by drivers that each act over a finite 
circular area, 1.91 cm (0.75 in.) in diameter, not a point. 

Microphone and Anemometer Operational Characteristics 

Figure 38 is a plot of the analytical and experimental values of 
time-mean anemometer output voltage E with respect to P^ ms > Tor various 
values of p 2 /P i - The experimental values of E are similar to those 
predicted analytical ly. 

Figure 39 is a plot of the analytical and experimental values 
of the root-mean-square pf the time-dependent anemometer output voltage 
*>ms respect to P^ ms for various values of P 2 /Pi- The analytical 

and experimental values are not similar. 

Figures 40, 41, 42 and 43 illustrate the analytical and experimental 
microphone and anemometer waveforms for the (100), the (202) and the 
combined (100) plus (202) spinning modes. The waveform for the combined 
mode depends on the phase relation <}> 2 . The waveforms for <j> 2 equal to 
zero and 180 degrees are shown. There are differences between the 
analytically predicted and experimentally obtained wave forms. The 
locations of the anemometer peaks with respect to the microphone peaks 
are not the same analytically and experimentally -- a phase shift 
exists. The relative amplitudes are also different. 

Figures 44, 45 and 46 are frequency spectrum analyses of the 
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Figure 40. Analytical and Experimental Microphone 
and Anemometer Waveforms for the (100) 
Spinning Mode 
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Experimental 


Figure 41. Analytical and Experimental Microphone 
and Anemometer Waveforms for the (202) 
Spinning Mode 
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Experimental 


Figure 42. Analytical and Experimental Microphone 
and Anemometer Waveforms for the 
Combined (100) plus (202) Modes with 
P 2 /P 1 = 0.6, 4 , i = 4 > 2 = ®i = ®2 = ^ 
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Figure 43. Analytical and Experimental Microphone 
and Anemometer Waveforms for the 
Combined (100) plus (202) Modes with 
P 2 /P 1 = 0.6, <h = = 62 = 0 and 4> 2 = 
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Frequency (kHz) 

Figure 44. Frequency Analysis of the Microphone and Anemometer Signals for the 
(100) Spinning Mode 



Relative Amplitude (dB) 



Frequency (kHz) 

Figure 46. Frequency Analysis of the Microphone and Anemometer Signals for the 
(100) plus (202) Combined Spinning Modes for p 2 /Pi = 0.6 
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waveforms for the (100), the (202) and the combined (100) plus (202) 
spinning modes. 

The (100) mode microphone signal is predicted to be a pure sine 
wave at f 100 . The frequency analysis of figure 44 indicates it has, 
in addition, a component at 2f 100 ; however, the amplitude of this 
harmonic component is about 25 dB below that of the f 100 so it is not very 
significant. The results of the frequency analysis of the (202) mode 
microphone signal, shown in figure 45, are similar, but the amplitude 
of the harmonic component is about 35 dB below that of the f 202 - 

The particular combined (100) plus (202) spinning mode shown 
in figure 46 had a ratio of p 2 /Pi of 0.6. For this ratio the sound 
pressure levels should differ by 4.44 dB, and this is approximately the 
value obtained from the frequency analysis shown in figure 46. 

The velocity variations associated with the (100) mode are 
predicted to be sinusoidal. The anemometer acts as a nonlinear 
rectifier. Rectification doubles the fundamental frequency. Hence, a 
frequency analysis of the anemometer signal should be composed of a 
strong component of 2f 100 plus higher harmonics. The frequency analysis 
of the anemometer signal presented in figure 44 shows that the expected 
frequencies do occur; but, in addition, there is a subharmonic component. 
Examination of the experimental anemometer waveform in figure 40 also 
indicates that its periodicity is associated with f 100 not 2f 100 . The 
same phenomenon occurs for the (202) mode of figure 45. 

The frequency associated with the anemometer signal of the 
combined mode should be f 100 plus higher harmonics, as is shown in 
figure 46. 
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Response Factors 

Figure 47 is a plot of the experimental values of the response 

factor as a function of P' „ for various amounts of harmonic content 

rms 

P 2 /Pi- A representative set of data at each value of p 2 /Pi is displayed. 

This figure very clearly shows the strong effect of the distortion 

component on the response factor. The overall shape of these curves is 

similar to that which would be expected if a through flow existed in 

the chamber rather than the no -forced -through -flow condition which 

existed (see figures 9 and 10). 

Figure 48 shows the effect on the response factor of changing 

the phase angle <t> 2 by 180 degrees. It is a graphical display of the 

experimentally determined response factors with respect to P‘ for 

rms 

p 2 /pj equal to 0.4 and 0.6. 

Figure 49 shows the variation in the response factor with respect 
to the phase between the fundamental (100) component and the distortion 
(202) component of the acoustic field for p 2 /p i = 0.6. This figure 
illustrates the strong effect of the phase of the distortion component 
on the response factor. Figure 50 is the analytically determined values 
of as a function of <t> 2 for the same conditions as existed in figure 
49. The value of <j> 2 at which the maximum response occurs is different 
experimentally than it is analytically. Also, the analytical magnitudes 
of the response are much larger than those obtained experimental ly. 





















Chapter 4 


DISCUSSION OF THE ANALTYICAL AND 
EXPERIMENTAL RESULTS 

INTRODUCTION 

This chapter is devoted to a discussion of the analytical and 
experimental results. The discrepancies between the two are discussed 
and an extension of the results to rocket stability is made. 

ANALYTICAL RESULTS 

This research began with an analytical study of the response of 
a constant-temperature hot-wire anemometer to sinusoidal and distorted 
sinusoidal acoustic oscillations. This analytical study showed that 
the response of the anemometer is similar in nature to the response of 
a vaporization-limited combustion process. In both processes the 
amount and phase of the distortion component, relative to the fundamental 
component, strongly affects the response of the system. For some 
conditions, the addition of distortion increases the response by more 
than an order of magnitude above that expected for the undistorted case. 
The response is an estimate, obtained from open-loop studies, of the 
closed-loop performance of the system. Unstable closed-loop performance 
is definitely predicted for some cases. 

The analytical study contains predictions of the anemometer's 
response both with and without a longitudinal component of steady flow. 
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The no-flow response is greater than the corresponding response with 
flow. This effect is especially strong at low sound pressure levels. 

At low values of P^ ms the response for the no-flow condition increases 
while the corresponding response with flow decreases and tends toward 
zero. At high values of P^ ms both the flow and no-flow situation tend 
to have the same behavior -- as P^ ms increases, the response decreases. 

EXPERIMENTAL RESULTS 

An experimental system was designed and built to produce for the 
hot wire the same acoustic environment that was studied analytically. 

The results of the studies made of the acoustic fields indicate that 
the desired acoustic fields are being generated, but that they include 
some distortion. Secondary effects due to the air's viscosity (acoustic 
streaming) are strongly indicated. 

With the experimental system, it is possible to determine a 
response factor for the anemometer when the hot wire is exposed to 
various acoustic fields. The results of the response studies show 
that both the amount and phase of the distortion component of the 
acoustic field, relative to the fundamental component, strongly affect 
the response of the system. For some conditions the addition of 
distortion increases the response by more than an order of magnitude 
above that obtained for the undistorted case. 

Discrepancies 

In spite of qualitative agreement, a number of major differences 
exist between theory and experiment; they are: 

® The variation in R^ as a function of P^ ms resembles that 
associated with the through-flow case, not the no-flow case. 
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© The magnitudes of R ^ are not as large as predicted. 

© E' is much smaller than predicted, 
rms 

@ The anemometer signal waveforms are not as predicted; there 
are both phase and amplitude distortions. 

© The maximum value of does not appear to occur at the value 
of <t> 2 predicted by analysis. 

As a group, the discrepancies suggest possible explanations of the 
results. 

Acoustic streaming . The shape of the versus P^ ms curves 
strongly suggests the presence of a constant flow over the wire. Hribar 
[8] has shown that when a standing (101) or (201) acoustic field exists in 
the air in an open ended cylindrical chamber, large standing vortices 
develop in the chamber. These vortices are not predicted by nonviscous 
theory but result when viscous effects are included in the governing 
equations. The vortices are a secondary flow effect. They do not 
significantly affect the time-dependent pressure and velocity distribu- 
tions predicted by nonviscous theory. The velocities associated with 
these vortices {acoustic streaming) are analytically predicted to be an 
order of magnitude less than that associated with the time-dependent 
velocities resulting from the nonviscous solution. There is also likely 
to be acoustic streaming associated with the presence of the hot-wire probe. 

The addition of a time-mean tangential acoustic streaming velocity 
to the velocity predicted by the nonviscous solution would have an 
effect similar to that associated with a flow through the chamber. The 
acoustic streaming velocity would add a time-mean tangential velocity 
to the velocity field; axial flow through the chamber adds a time-mean 
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longitudinal component to the velocity field. The streaming velocity 

would alter the nature of the response so that R „ as a function of P* 

nl rms 

would take on the general shape associated with through-flow velocities. 

Temperature effects . If a standing (100) acoustic field were to 
be set up in the chamber, and a hot-wire probe placed at a velocity node, 
the anemometer should not detect a time-dependent velocity. When this 
is actually carried out, an apparent low level velocity is detected at 
the (100) frequency. As stated earlier, velocities are rectified by the 
hot wire; an oscillating velocity associated with the (100) mode is 
expected to have a fundamental frequency of 2f 100 . 

For a standing (100) acoustic field, a pressure, and therefore 
a temperature antinode exists at a velocity node. As the anemometer 
equations indicate, the anemometer is sensitive to temperature fluctua- 
tions. Because of this, what may look like a velocity fluctuation may 
actually be a temperature fluctuation. A temperature fluctuation due 
to a (100) acoustic field would occur at f 100 . 

If the hot wire were actually sensing the time-dependent tangential 
velocity associated with the spinning (100) acoustic field, plus a 
time-mean tangential velocity due to the acoustic streaming, plus an 
oscillating "velocity" due to temperature fluctuations, and the relative 
magnitudes of these were of the proper value, the output of the anemometer 
would have a waveform similar to that actually obtained for the (100) 
spinning acoustic field (see figure 40). 

Effect of E on Ep ms - The value of E^ ms is affected when a constant 
flow exists over the wire because of the nonlinear nature of the 
anemometer. For example, referring to figure 34, the change in E 
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resulting from a change in velocity from zero to one meter per second is 
greater than the change in E resulting from a change in velocity from 
one to two meters per second. As the constant velocity increases, the 
value of E^, ms associated with a fluctuating velocity, whose magnitude 
is remaining constant, decreases. 

Equipment limitations . The operating characteristics of the 
equipment also affects the microphone and anemometer outputs and, 
therefore, the response. Phase shifts in the microphone could cause 
the signal at one frequency to be shifted relative to the signal at 
the other; hence, the microphone signal which appears on an oscillo- 
scope to be that associated with <j> 2 equal to zero may actually be 
associated with another value of $ 2 - This effect could be why the 
maximum value of the experimentally determined response occurs at 
<t> 2 . ~ ~ 45° (figure 49) when analytically the maximum value of the 
response occurs at <j> 2 = 0 (figure 50). Some experiments on the micro- 
phone system were carried out to determine the phase shifts associated 
with the frequencies of interest; however, nothing conclusive was 
obtained from these investigations. 

The response of the anemometer may also be inadequate at the low 
levels of velocity being considered. The anemometer instruction manual 
[21] contains a graph which shows the amplitude and phase lag as a 
function of frequency for various time-mean velocities. Using this 
graph, it can be estimated that there is a phase lag of about ten 
degrees for a frequency of 1000 Hz and a phase lag of about twenty 
degrees for a frequency of 2000 Hz. When these phase lags are converted 
to time delays, the delay time for each signal is the same; hence, the 
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signal is delayed in time, but not distorted. This time delay effect 
would cause the peaks of the anemometer signal to be shifted relative 
to the peaks of the microphone signal. The same graph indicates there 
may be some amplitude distortion; however, the lowest velocity listed 
is ten meters per second so an actual value cannot be estimated with 
accuracy. 

The combined effect of both microphone and anemometer phase 
shifts and/or time lags would result in an uncertainty about the 
value of <J> 2 , the possibility that the microphone and anemometer signals 
being correlated do not originate at the same time and distortion of the 
P' and V waveforms. All of these factors would tend to reduce the 
resulting value of R 

Limitations of the theory . In addition to all the above factors, 
it is also a possibility that the theoretical development itself is not 
entirely valid. In theory the velocity passes through zero twice during 
each cycle. The "zero" velocity region during steady-state, steady- 
flow operation is free convection dominated. The equations were developed 
for forced \convect ion heat transfer from the wire. Also, the assumption 
was made that the steady-state, steady-flow calibration equation of a 
hot wire could be applied to an unsteady situation if the instantaneous 
values of the fluid properties were used on the calibration equation. 

This assumption may not be valid for this extreme application of the 
hot-wire anemometer. 

The above discussions present possible explanations of the 
experimental observations. They are aspects of this study which should 
be explored during future investigations. 
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S ignificance to Rocket Stability 

The results obtained indicate why, when an instability occurs 
in a rocket engine, a sinusoidally varying time-dependent velocity 
component of the sound-pressure field very rapidly becomes characterized 
by a fundamental frequency distorted by the second harmonic frequency. 

The frequency spectrum analyses of figures 48, 49 and 50 indicate 
how a feed back mechanism may develop which will allow an initially 
sinusoidal disturbance to become a distorted sine which can increase in 
intensity very rapidly. First, the chamber will respond much more 
strongly to its resonant frequencies than to other frequencies. (This 
was demonstrated by the resonant frequency analysis performed as part of 
this research.) It is able to pick out and respond to its resonant 
frequencies from a spectrum of frequencies which may exist in the 
chamber. When the chamber responds to fi 00 , a signal with a strong 
2f 100 frequency component results from the anemometer. If this signal 
is fed back into the chamber, the chamber will respond strongly to it 
since it is at the (202) resonant frequency. The anemometer output 
resulting from the combined (100) plus (202) acoustic fields has strong 
components at both the f 100 and 2f 100 frequencies. If this signal is 
then fed back into the chamber, the intensity of both the (100) and the 
(202) acoustic fields will be amplified. This loop, under proper 
feed-back conditions, would very quickly build up the sound field to a 
very high intensity. 

Whatever secondary effects are found to occur in the chamber studied 
here would also tend to occur in a rocket combustion chamber. The 
proposed description of the feedback mechanism for amplifying an instability 
would also apply to the closed-loop situation in a rocket engine. 
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CONCLUSIONS AND RECOMMENDATIONS 
INTRODUCTION 

An analytical and experimental investigation of the effect of 
sinusoidal and distorted sinusoidal acoustic perturbations on a Reynolds 
number dependent process has been made. In both the analysis and the 
experiment, the Reynolds number dependent process examined was the 
response of a constant-temperature hot-wire anemometer. Before the 
experimental investigation of the anemometer's response could be 
carried out, an experimental system had to be developed and tested to 
determine if it could produce for the hot wire the same environment 
examined analytically. 

Since there were differences between the analytical and experi- 
mental data, some recommendations are made for further studies which may 
clarify the reasons for the differences. In addition, recommendations are 
made for extensions of this research. 

CONCLUSIONS 

1. The results of the analytical investigation of the response 
of a constant-temperature hot-wire anemometer to a sound field show 
that the anemometer's response is qualitatively the same as that predicted 
for vaporization limited combustion -- the addition of distortion to 
sinusoidal perturbations of the environment of a hot wire can increase 
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the response significantly, and the relative magnitude and phase of the 
distortion component are predicted to be the dominant variables in 
determining the response. 

2. The experimental investigation shows that it is possible to 
set up a physical system which can generate the same acoustic environment 
for a hot-wire probe that was investigated analytically and which can 

be used to determine the response of the hot-wire bridge voltage to 
that environment. 

3. Investigations performed with the experimental system show 
that the relative amplitude and phase of the distortion component of 
the hot wire's environment very strongly affect the response; the 
distortion component can increase the response by an order of magnitude 
above that for the undistorted case. 

4. The distorted acoustic fields generated cannot be totally 
characterized by the equations resulting from the solution of the 
inviscid wave equation. Viscous effects should be included in the 
analysis in order to account for acoustic streaming and any other 
secondary phenomena which may occur. 

5. The assumption that the steady-state, steady-flow calibration 
equation of a hot wire can be applied to an unsteady situation, if the 
the instantaneous values of the fluid properties are used in the 
calibration equation, may not be valid. 

6. A feedback and amplification process is possible when two 
resonant frequencies exist in a chamber such that the one resonant 
frequency is twice that of the other. 
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RECOMMENDATIONS 

1. More extensive investigations of the acoustic fields in a 
cylindrical chamber need to be performed with special attention being 
given to acoustic streaming phenomena and temperature fluctuations. 

2. An investigation of the limitations on the assumption that the 
steady-state, steady-flow calibration equation of a hot wire can be applied 
to an unsteady situation, if the instantaneous values of the fluid properties 
are used in the calibration equation, needs to be carried out. 

3. The dynamic behavior of the anemometer needs to be examined in 
more detail, especially with regard to E^ ms - 

4. The dynamic behavior of the microphone needs to be examined 
in more detail. 

5. More data of the type reported herein need to be taken and 
analyzed, 

6. Data should be taken for other wire locations and orientations 
to establish how the response varies with chamber location. 

7. An analytical study of the effect of various length to 
diameter ratios on the resonant frequencies of the chamber should be 
performed since the resonant frequencies that are multiples of lower 
mode frequencies affect the response. 

8. A theoretical analysis needs to be done which will Quantitatively 
relate the experimental results obtained for the constant-temperature 
hot-wire anemometer to the results to be expected for vaporization- 
limited combustion. 

9. A closed-loop system should be designed with which it can be 
determined if the open-loop predictions of the closed-loop performance 
are valid. 



REFERENCES 


140 



LIST OF REFERENCES 


1. Priem, R. J. and M. F. Heidmann, "Propellant Vaporization as a 

Design Criterion for Rocket Engine Combustion Chambers," NASA 
TR R-67 , 1960. 

2. Heidmann, M. F. , "Amplification by Wave Distortion in Unstable 

Combustors," AIAA Journal , Vol. 9, No. 2, Feb. 1971, p. 336. 

3. Heidmann, M. F. , "Amplification by Wave Distortion of the Dynamic 

Response of Vaporization Limited Combustion," NASA TN D-6287 , 
May, 1971. 

4. Sotter, J. C. and G. A. Flandro, "Resonant Combustion in Rockets," 

Sci. Am ., Vol. 219, No. 6, Dec. 1968. 

5. Conrad, E. W. , H. E. Blommer, J. P. Wanhainen, and D. W. Vincent, 

"Interim Summary of Liquid Rocket Acoustic-Mode-Instability 
Studies at a Nominal Thrust of 20,000 Pounds," N ASA TN D-4968 , 
Dec. 1968. 

6. Harrje, D. T. and F. H. Reardon, editors, Liquid Propellant Rocket 

Combustion Instability , NASA SP-194, 1972. 

7. Rayleigh, J. W. S., The Theory of Sound, Dover, New York, 1945, 

Vol. 2, p. 226. 

8. Hribar, A. E. , "The Effect of a Transverse Resonant Acoustic Field 

Upon Viscous Fluid In a Cylindrical Enclosure," Ph.D. Thesis, 
Purdue University, 1969. 

9. Osborn, J. R. and J. M. Bonnell, "Transverse Modes of Acoustic 

Oscillations," Jet Propulsion Center, Purdue University, Report 
No. TM-64-4, June 1964. 

10. Osborn, J. R. and J. R. Rahon, "Effects of Radial Energy Releast 

Variations on Transverse Combustion Pressure Oscillations," 

Jet Propulsion Center, Purdue University, 1-62-3, March 1962. 

11. Heidmann, M. F. and P. R. Wieber, "Analysis of Frequency Response 

Characteristics of Propellant Vaporization," NASA TN D-3749, 
Dec. 1966. 


141 



142 


12. Priem, R. J. and D. C. Guentert, "Combustion Instability Limits 

Determined by a Nonlinear Theory and a One-Dimensional Model," 
NASA TN D-1409 , 1962. 

13. Priem, R. J. , "Influence of Combustion Process on Stability," 

NASA TN D-2957 , Aug. 1965. 

14. Purdy, K. R. , M. B. Ventrice and J. C. Fang, "An Analytical 

Investigation of the Open-Loop Amplification of Reynolds Number 
Dependent Processes by Wave Distortion," Proceedings of the 
Tenth Annual Southeastern Seminar on Thermal Sciences , 1974. 

15. King, L. V., "On the Convective Heat Transfer from Small Cylinders 

in a Stream of Fluid. Determination of Convective Constants 
of Small Platinum Wires with Application to Hot-Wire 
Anemometry," Phil. Trans. R. Soc ., Vol . 21 4A, 1914. 

16. Andrews, G. E. , D. Bradley and G. F. Hundy, "Hot-Wire Anemometer 

Calibration for Measurements of Small Gas Velocities," Int . 

J. Heat Mass Transfer , Vol. 15, 1972, pp. 1765-1785. 

17. Collis, D. C. and M. J. Williams, "Two Dimensional Convection from 

Heated Wires at Low Reynolds Numbers," J . of Heat Transfer, 

Vol. 6, 1969, pp. 357-384. 

18. Hilsenrath, J., et. al., Tables of Thermal Properties of Gases, 

National Bureau of Standards, Circular 564, Nov. 1955. 

19. Purdy, K. R. , G. W. Simmons, A. E. Hribar and E. I. Griggs, 

"Acoustically Augmented Convective Drying," Sonochemical 
Engineering , Vol. 67, 1971, pp. 55-67. 

20. Fang, J. C.and K. R. Purdy, "Calibration of a Constant-Temperature 

Hot-Wire Anemometer for Very Low Reynolds Numbers," ME-44 Report, 
Department of Mechanical Engineering, Tennessee Technological 
University, 1973. 

21. DISA, "Constant Temperature Anemometer Instruction Manual for Type 

55A01 . " 

22. Purdy, K. R. , "Viscous Fluid Flow Under the Influence of a Resonant 

Acoustic Field," Ph.D. Thesis, School of Mechanical Engineering, 
Georgia Institute of Technology, 1963. 

23. Morse, P. M. , Vibrations and Sound , Second edition, McGraw-Hill 

Book Co., Inc., New York, 1948. 

24. Abramowitz, M and I. A. Stegun, editors, Handbook of Mathematical 

Functions , National Bureau of Standards, Applied Mathematics 
Series No. 55, 1964. 

25. DISA, "Instruction and Service Manual for Type 55D41/42 Calibration 

Equipment," 



APPENDICES 


143 



APPENDIX A 


SOLUTION OF THE INVISCID WAVE EQUATION 
PROBLEM STATEMENT 

The system being analyzed is shown in figure A1 . 


z = 0 z = L 



Figure A1 . System Model 

It consists of a cylinder of length L and radius R which is closed at 
both ends and contains an ideal gas. The following assumptions are 
made about the system: 

© The gaseous medium is in a uniform state except for 
weak isentropic fluctuations in the state properties 
and velocity caused by the acoustic field. 

© There is no heat transfer. 

© There is no mass flow into or out of the cylindrical 
enclosure. 
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The governing equations are: 
© Continuity 
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5p 


+ 


7-(p V) = 0 


(A. 1 ) 


® Momentum 


pj^=-?P (A.2) 

® Equation of State 

P - p R g T {A. 3} 

© Second Law of Thermodynamics 

s = constant (A. 4) 

Since the gas is inviscid, a velocity potential can be defined 

v$> = V (A. 5) 

Introducing this definition into the governing equations yields a non- 
homogeneous wave equation for the velocity potential 

V 2 * - ^2 f^T = g(*) (A. 6) 

oo 

It can be shown that the terms composing g(<j>) are at least an order of 
magnitude smaller than the left-hand side of equation (A. 6). A detailed 
order of magnitude analysis demonstrating this for the one-dimensional case 
is given by Purdy [22] and a further discussion can be found in Morse [23]. 

In accordance with the first assumption, the velocity, density 
and pressure can be expressed as the sum of a time-mean and a time- 
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dependent component as follows: 

V(r,0,z,t) = 7(r,e,z) + V'(r,e,z,t) 

P (r,e,z,t) = p"(r,e ,z) + p'(r,e,z,t) 

P(r,e,z,t) = P”(r,e,z) + P'(r,e,z,t) (A. 7) 

where the bar overscore denotes the time-mean value and the prime 
denotes the deviation from the mean value. 

The isentropic speed of sound for an ideal gas is 

= Y R g T = Y P/p (A. 8) 

The speed of sound in equation (A. 6), c ot , is based on the "undisturbed" 
property values -- the property values when V(r,e,z,t) = 0. 

Thus 



c 2 = k P /p 

OO CO oo 

(A. 9) 

The velocity potential is related to 

the acoustic velocity, pressure 

and density such that 




V' = va> 

(A. 10) 

n. - a* 

p • • p - n 

(A. 11) 

i _ P « 8$ 

p " c 2 at 

CO 

(A. 12) 


The boundary conditions are: 

© Acoustic pressure P 1 is finite at r = 0, or 
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),6,z,- 

at 


finite 


(A. 13) 


• Radial acoustic particle velocity is zero at the cylinder 
walls (r = R), or 


3»(R,e,z,t) . « 
ar 


(A. 14) 


For closed ends the longitudinal acoustic particle velocity 
is zero at the cylinder ends (z = 0,L), or 


a$(r,8,0,t) _ a$(r,e,L,t) _ 0 
az az 


(A.15) 


In cylindrical coordinates, the wave equation to be solved is 


1 3 ^ r ar^ + „1 a 2 $ , 1_ a£$ _ n 

r ar r 2 az 2 " c 2 at 2 


(A. 16) 


Assuming that the velocity potential can be expressed as 


<j>(r,e ,z,t) * Real [<j>(r,a ,z)e lwt ] 


(A. 17) 


substituting equation {A. 17) into equation (A. 16) and simplifying, yields 
the Helmholtz equation for <ji(r,0,z). 


v 2 <j> + k 2 <j> - 0 


(A. 18) 


where 


k 2 = u 2/c 2 


(A. 19) 


Equation (A. 18) can be solved using the method of separation 
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of variables. Assuming that 

4 >(r,e,z) = r(r,e)A(z) (A. 20) 

substituting equation (A. 20) into equation (A. 18), separating variables 
and using the arbitrary constant k z 2 , yields the following two equations: 

0 + k 2 *A = O (A.21) 

( r 1L) 

r + (k 2 - k z 2 )r 2 r = 0 (A. 22) 

The solution of equation (A.21) is 

a{z) = A cos(k 2 z) + B sin(k 2 z) (A. 23) 

Equation (A. 22) is a function of two variables. Letting 

r(r,e) = n(r)T(e) (A. 24) 

Substituting equation (A. 26) into equation (A. 24), separating variables 
and using the arbitrary constant m 2 , yields 

+ m 2 T = 0 (A. 25) 

and 

fpr+flf + < k2 - k z 2 - ■ 0 (A.26) 

The solution of equation (A.25) is 

T(e) = C cos(me) + D sin{me) (A. 27) 
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where m can take on only integer values since otherwise t(9) would be 
multiple valued. 

Equation (A. 26) can be put in the form of Bessel's equation. 

Letting 

k r 2 » k 2 - k 2 2 
r » x/k r 

y « Q{r) (A. 28) 

equation (A. 26) becomes 

x 2 |^r + x + (x 2 - m 2 )y = 0 (A. 29) 

which is Bessel's equation of order m. The solution is [24]: 

y = F J B (x) + GY r ( x) for m an integer (A. 30) 

where 


J m (x) = Bessel functions of the first kind of order 
m and 

Y m (x) = Bessel functions of the second kind of order 
m or Weber functions 

Converting back to the original variables, equation (A. 30) becomes 

fi(^) = F J m (k r r) + G Y m (k r r) (A. 31) 

Combining equations (A. 17), (A. 20), (A. 23), (A. 27) and (A. 31) 


gives 
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$(r,o,z,t) = Real [F J m (k r r) + 6 Y m (k r r)][C cos(mo 

+ D sin(me)] * [A cos(k z z) + B sin(k z z)]e iwt j 

(A- 32) 

Application of the first boundary condition, equation (A.15), 

requires that G = 0 since Yfk r) -*• - ® as r + 0. 

rn r 

Application of the second boundary condition, equation (A. 16), 
requires that 

J m '(k r R) = 0 (A. 33) 

J m ' is almost a periodic function; therefore the eigen values of equation 
(A. 33) will be defined as a mn u, 

k R = a it 
r mn 

k r = a mn 7 t/R ; m,n = 0,1,2,... (A. 34) 

For each value of m, there is a series of solutions which satisfy equation 
(A. 34); the symbol n is used to distinguish these. 

Application of the third boundary condition, a pressure node at 
z = 0,L, requires B = 0 and sin(k z L) = 0 — a solution can be obtained 
only for the eigenvalues 

k z = n z tt/L ; n z = 0,1,2,... (A. 35) 

Substituting the above results into equation (A.32) and combining 


constants yields 
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4 >(r,e,z,t) = Real (j m (a mn ,r/R) C cos(me) + D mnn sin(me)]- 


cos(n irz/L)e^ wt } 


(A. 36) 


Extracting the real part of equation (A. 36) and recognizing that for a 
linear equation the sum of the solutions is also a solution and that 
depends on the values of m, n and n » the general solution of the 
inviscid wave equation is 


*(r,e,z,t) = z J m^mn^ r/R ^ C0S ^ n z 1rz/L ^ A mnn cos ( me 
m,n,n, z 


+ (a t+r ) + B cos(me-oj t-n )] (A. 37) 

mnn 2 ^mnn z mnn z ' mnn z 'mnn /J v ‘ 


where A m „„ and 13 „ are constants which are associated with particular 

mnn^ mnn^ r 

values of m, n and n„ and r, „ and n „ are arbitrary phase angles also 

z mnn z mnn z 

associated with particular values of m, n and n .• 

The resonant frequencies of oscillation are specified by the 
eigenvalues of the solution. Using the relation to = 2irf, equations 
(A. 19) , (A. 28), (A. 34) and (A. 35) yield 



(A. 38) 


This specifies all possible resonant frequencies. 

The solution, equation (A. 37), can be put in the form of a standing 


wave or a spinning wave. 
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STANDING WAVE SOLUTION 


In equation (A. 37} the last term on the right-hand side is the 
sum of two traveling waves, one moving in the positive e direction, the 
other moving in the negative 0 direction. To form a standing wave the 


amplitudes of the two must either be equal, A = B with c 

mnn z mnn z mnn 2 

or of opposite siqn, A = B with c - n + ir. 

mnn r ~ 73 7 mnn mr \in - *' 


mnn z mnn z 


mnn z mnn z ' 


Assuming A mnn = B mnnz , equation (A. 37) can be put into simpler form, 
namely 1 


t(r,e,z,t) = i G J m (a ,r/R)cos(n ,z/L) 
mnn_ z 


cos(me)cos{w t + p ) 
' mnn z mnn z ' 


where G mnn are cons tants and the are arbitrary phase angles 

associated with particular values of m, n and n z< 

Using equations (A. 39) and (A. 11) a general term for P’ is 


P ™nn z = s »“rann z G rann z J m ( ^WRjcosfn^z/L) 


cos(me)sin(ui ^ t + g ) 
mnn z mnn z 


Letting the maximum value of P^ nn for a particular {mnn ) mode be 

^ z 

P m „„ , i.e., 
mnn z 


^mnn z p «> w mnn z G mnn z ^m^mn 1 ^ 


solving equation (A. 41) for G 
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mnn 

-v _ Z 

3 mnn p uT~ J ( a tt } 
z °» mnn, m mn 


(A. 42) 


and substituting equation (A. 42) into equation (A. 39), but considering 
only the general term, gives 


$ n (r,9,z,t) = 
mnn, p 

z <*> mnn_ m mn 


mnn 7 J m^ a mn irr/R ^ 

2 rn HI — , — cos(n 7 rz/L) 
a it) Z 


cos (me) cos (ui, t + e ) 
mnn z mnn^ 


(A. 43) 


Using equation (A. 43), general terms for the acoustic pressure, velocity 
and density are: 

« J_(a__7rr/R) 

P'n = p n -i / - — v — cos(n_ifz/L)cos(me)sin(ai t 
mnn z mnn z Wl z rann z 


+ $ ) 
mrm ' 


mnn. 


(A. 44a) 


V r 


mnn, a mn ¥ J* (a ttt/R) 


nr mn 


0 it) 

mnn z ninn z 


r J m( a mn 7r ) 


cos(n irz/L)cos(me} 


m v mn 


cos(w t + 3 ) 

' mnn z mnn^ 


(A. 44b) 


r e. 


P mnn, m J (cc Ttr/R) 
z nrmn 


mnn 


P w , r J m (a m „Ti) 
00 mnn m mn 


cos(n z 7rz/L)sin(me) 


COS (a> t + 8 „ ) 

mnn z mnn z 


(A.44c) 
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P f nnn n z * J (a mn irr/R) 

’ = i , / mn — < — sin(n 7rz/L)cos(me) 

z_,_„ p L J„ a z ' ' 

mnn » mnn m mn 


c ° s ( (a) t + e „ ) 

mnn z mnn^ 



SPINNING WAVE SOLUTION 

As stated earlier, the general solution of the inviscid wave 
equation, equation (A. 37), is actually the sum of two traveling waves -- 
one moving in the positive e direction, the other in the negative e 
direction. Since the waves move tangentially the terminology "spinning" 
is used. Setting equal to zero in equation (A. 37) yields a 

solution composed of traveling waves which spin in the positive 
direction, namely 


$ (r,e,z,t) = £ J m (a mn iTr/R)cos(n 7 TTz/L) 


cos (me - io t - n ) 
mnn z mnn z 'mnn z ' 


In this form the interchangabil ity of the effect of changes in either 
time or angular position can be seen. 

Trigonometric transformation of the last term of equation (A. 45) 


yields 
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$(r,e ,z,t) z J (a Trr/R)cos(n z 7rz/L) 
mnn 2 


B [cos(me)cos(oj „ t + n„ _ 

mnn„ mnn_ mnn.. 

Z Z i 


- sin(me)sin(a) 


mnn 


t + 


n„„„ )] 

mnn z 


(A. 46) 


Equation (A. 46) is in the form of the sum of two standing waves; hence, 
a spinning wave can be considered to be composed of two standing waves 
of equal amplitude, but 90/m degrees out of phase with each other in 
angular position and 90 degrees out of phase with respect to time. 

This is a useful concept for experimental work. 

Following a procedure similar to that used for the standing case, 
the general term for the velocity potential is 


4> 

mnn z 


(r,e,z,t) 


z m x mn ' 


p w _ T ( oTT tT) 

oo mnn m mn 


cos(n z 7rz/L) 


cos (me - u> mnn t - n „ ) 
mnn z mnn z 


(A. 47) 


and the general terms for the acoustic pressure, velocity and density 
are: 


- J (a nr/R) 

P' . = P mnn — t — t y — cos(n irz/L)sin(me 

mnn. mnn. J (a tt l z 

z z m mn 


- a) t n ) 
mnn z 'mnn^ 


(A. 48a) 
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P mnn J ' (a Trr/R) 

V' = p i j — • — \ cos(n irz/L)cos(me 

r p R J tci tt) 2 

mn z °° mnn z m mn 


- (jj t - n ) 

mnn z mr 1 ' 


v ; 


mnn. 


P mnn, m J m (a m „WR) 

— — tj—j y— cos(n irz/L)sin(me 

p ~ w mnn_ r VW* z 


— aj t - n ) 

mnn z Win ' 




^mnn, 1 ^ J (a irr/R) 
z m mn 


mnn. 


P cd L J m (a m „ir ) 
« „„„ m mn 
mnn. 


sin(n z Tiz/L)cos(m0 


■a) t ” n ) 

mnn z 'mnn 7 


mnn, J {01 -rrr/R) 
z m mn ' 


mnn c 2 J (a tt) 

z °° m' mn 


cos(n z TTz/L)sin(m0 


- 03 t - n ) 

mnn z , mnn z / 


(A. 48b) 


(A. 48c) 


(A.48d) 


(A.48e) 



APPENDIX B 


SELECTED VALUES OF « 

mn 

The resonant frequencies of a closed-ended cylindrical chamber 
depend on the eigenvalues of equation (A. 33), 


J m '(k r R) - 


(A. 33) 


Since J * is almost a periodic function, the eigenvalues of equation 
(A. 33) are defined as a^ir; thus 


k 


r 



7T 


Substituting equation (A. 34) into equation (A. 33) gives 


(A. 34) 


J ' (« TT ) = 0 

m mn 


(B.l) 


For each value of m there is a series of solutions which satisfy equation 
(B.l) — the symbol n is used to distinguish these. Table B.l lists 



c\j co lo <£> oo 
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Table B.l. 


m n=0 n = 1 n = 2 n = 3 n = 4 


0 0.0000 1.2197 2.2331 3.2383 4.2411 

1 0.5861 1.6970 2.7140 3.7261 4.7312 

0.9722 2.1346 3.1734 4.1923 5.2036 

1.3373 2.5513 3.6115 4.6428 5.6624 

1.6926 2.9547 4.0368 5.0815 6.1103 

2.0421 3.3486 4.4523 5.5108 6.5494 

2.3877 3.7353 4.8600 5.9325 6.9811 

2.7304 4.1165 5.2615 6.3477 7.4065 

3.0709 4.4931 5.6576 6.7574 7.8264 




APPENDIX C 


HOT WIRE CALIBRATION RESULTS 


The equation which relates the velocity of flow over a hot wire 
to the voltage output from a constant temperature hot-wire anemometer 
associated with it is given by equation (2.17) 


r 


E = 


[A + B Re- N (V 2 ) N/2 ] 
'1 


1 

0 . 17 


T 


T 

L 00 J 


. °J 


1 0*877 


[T u - T] 


} 


0 ' 5 


(2.17) 


For the steady-state, steady-flow conditions which exist during wire 
calibration, all quantities in equation (2.17) are constant except the 
velocity; hence, equation (2.17) reduces to 

~M °* 5 

E = [A + E V N ] (2.18) 


where 


and 


A - A 


[V 

0.17 

T f 

T 


T. 

coj 


L o j 


0.877 


£ T w - 


(2.19a) 


B = B Re- 
T 1 


N 

0.17 

— f 
-h 

t 

T 


T 

oo 

„ - 


oj 


0.677 


[T w - T] 


(2.19b) 


For steady-state, steady-flow conditions, T = T and P = P . 

oo oo 

During the calibration of a wire, the voltage output from the 
anemometer and the corresponding pressure drop across the nozzle of the 
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wind tunnel in which the wire is located is obtained. The equipment 
and procedure used to do this is described in Chapter 3. For the 
equipment described, the pressure drop across the nozzle of the wind 
tunnel is measured as the amplified analog voltage of an MKS Baratron 
pressure transducer. The voltage so obtained is converted, with a 
calibration constant determined by experiment, to a pressure difference. 
Table C.l lists the data and resulting pressure difference obtained 
during a calibration of the wire. The calibration constant which 
relates H to aP is 


H/aP = 122.53 vol ts/inch-of-water 
Table C.l. Hot-Wire Calibration Data 


Anemometer Output 
E 

(volts) 

Baratron Output 
H 

(volts) 

Pressure Drop 
AP 

(in. of H 2 0) 

Velocity 

V 

(m/sec) 

3.8160 

0.0314 

0 

0 

4.0235 

0.0472 

0.000129 

0.235328 

4.1374 

0.0722 

0.000334 

0.378676 

4.2289 

0.1048 

0.000601 

0.507979 

4.2973 

0.1386 

0.000878 

0.613984 

4.3823 

0.2055 

0.001427 

0.782736 

4.4535 

0.2854 

0.002082 

0.945472 

4.5156 

0.3825 

0.002878 

1.111608 

4.5755 

0.4992 

0.003834 

1.283022 

4.6397 

0.6392 

0.004982 

1.462545 

4.6975 

0.8254 

0.006508 

1.671596 

4.7459 

0.9882 

0.007842 

1 .834936 

4.7864 

1.1914 

0.009508 

2.020472 

4.81 77 

1.2817 

0.010248 

2.097624 


The velocity of the air passing over the wire in the nozzle of the 


wind tunnel can be determined with the following formula [25]: 
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2y aP (Y ' 1)/Y 

V= (^T )R aTH-O-^ 3 (3.13) 

This velocity is also shown in Table C.l. 

Using the information given in Table C.l and the least-squares 
curve-fitting procedure described in Chapter 3, the constants A, B and 
N are found to be 

A = 11.909 
B = 9.0995 

N = 0.444 

with a goodness-of-fit of 0.99988. (The data at a velocity of zero are 
not used since the equations were developed for forced-convective 
heat transfer and, at zero velocity, free-convective heat transfer 
is dominant.) 

The following additional data were also recorded: 

R » 4.02 n 

oo 

R w = 5.42 a 

T = 275°K 

00 

P = 736.79 mm Hg 

oo 

D « 5 x 10" 6 m 
w 

The temperature of the wire T w can be found from 

+ (2.7) 

where a = 0.0126 fl/K°. 

The film temperature is determined from its definition 



T f = < T w * T)/2 

The quantity Re* is 

T i 


Re, 


p Vi D, 


w 


R a T f 


0.757 


From the ideal gas equation 
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(C.l) 


(2.14) 


p = (2.10) 

where 

R = 287.1 m 2 /sec 2 -K° 

a 

From Hilsenrath [18] 

u 0 = 1.716 x 10" 5 Kg/m-sec 

Combining equations (2.19) with equations (2.7), (C.l), (2.14), (2.10) 
and the values of A, B, R , R , T , P , D, a, R, and the values of 

oo yf oo co yj a Q 

A A 

A and B become 

A = 0.5205 
B = 0.6729 


Substituting A, B and N into equation (2.17) yields 

r -.o - 17 


- j[0. ! 


lr 0.444 ~ 0.222, 

E = <[0.525 + 0.6729 Re f (V 2 ) ] 


Ft, 


L T oJ 


0. 877 


0.5 


[T w - T] 


(C.2) 
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